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The characf eristic equation for the propagation constants of waves in a 
filled circular guide of arbitrary radius is written in terms of magnetizing 
field and a carrier density, which are shown essentially to determine the 
dielectric and permeability tensors for a gas discharge plasma and for a 
feriiie. The complex structure of the spectrum, of propagation constants and 
its dependence upon radius and the two parameters are analyzed by a semi- 
graphical method, supplemented by exact formulae in special regions. Thus 
the course of individual modes may be charted with fair accuracy, 

1. INTRODUCTION 

Any material medium which propagates electromagnetic disturbances 
possesses a kwal electric or magnetic stmcture and it is just the motion of 
the electric or magnetic carriers imder the fields of the disturbance that 
determines how the propagation takes place. If a dc magnetic field be 
appHed to the mediiun one may expect the local response to be altered 
and, consequently, to find changes in the character of the propagation. 
Oyrumagnetic media arc those for wliich such changes are suflficiontly 
large to be experimentally significant. For plane waves and for optical 
frequencies the experimental effects and their explanation have been 
familiar for a great many years. The non-reciprocal rotation of the plane 
of polarization of light travelling parallel or antiparallel to an applied dc 
magnetic field, which is known as the Faraday effect, is such a phenom- 
enon. So also is the fact that the medium becomes doubly refracting for 
arbitrary directions of propagation. 

Interest in gyromagnetic media at longer wavelengths first arose in 
connection with radio propagation in the ionosphere. The ionosphere is 
essentially an ionic cloud and the earth supplies a magnetic field, which, 
for the charge densities involved, is sufficient to produce a large effect 
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upon propagation. Here, as in the earlier optical cases, the disturbances 
considered are essentially plane waves. In recent years, with the ex- 
tensive development of microwave techniques, two gyromagiietic media 
have been investigated using guided waves. One of these is the gas dis- 
charge plasma, an ionic medium like the ionosphere, in which, however, 
the charge density may be varied over wide ranges in a controllable 
manner. The magnitude of the effects observed in such ionic media are 
governed by the relation of the apphed frequency to the cyclotron fre- 
quency of the ions in the dc magnetic field. Goldstein and his associates' 
have studied the propagation of waves in a cylindrical waveguide within 
which a discharge is supported and to which a longitudinal magnetic 
field is applied. iVmong many effects which they have observed is a 
large Faraday rotation. 

The other medium being actively investigated is the low-loss ferro- 
magnetic medium, as exemplified by the ferrites. In this case the pe- 
culiarities of the medium have their origin in the precession of the 
magnetization of the ferrite about the applied field. This precession 
takes place with a frequency dependent upon the applied field strength 
and large changes in the nature of the propagation occur when the fre- 
quency of the r.f. applied field approaches this. Polder- woi'ked out the 
effective properties of such a medium for plane waves and Hogan''' has 
made various experimental studies of the propagation in cylindrical 
guides containing ferrite. Here, again, Faraday rotation and other 
non-reciprocal effects have been observed. 

In this paper a variety of topics associated with the theory of guided 
waves in gyromagnetic media is considered, with the main emphasis laid 
on the ferrites. The exposition does not attempt to be systematic. Very 
few problems in this field admit of a thorough analytic treatment and, 
frequently, the more closely allied they are to the practical uses of ferrites 
in microwave devices the more fragmentary is the analysis. On the other 
hand since the problems can always be formulated it is always possible 
in specific cases to resort to a purely numerical solution. The problems 
considered here all arise in the effort to analyze the operation of \'arious 
devices and different idealizations are utiHzed in particular cases. 

In Part I the general properties of gyromagnetic media are discussed 
and the connection between the phcnomenological constants of the 
medium and the underlying molecular model is derived for the ferrite 
and for the plasma. The assumptions necessary to render the ferrite 
problem tractable are discussed at some length. Maxwell's equations are 
written down for a general gyromagnetic medium and some of the salient 
features of their solution are noted. The propagation of circularly po- 
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larized waves in circularly cylindrical guide filled with ferrite or plasma 
is then considered. The characteristic equation connecting frequency 
and propagation constant is first derived. For the purpose of obtaining 
results which can be compared with experiment, a specific molecular 
model is chosen for the ferrite. In this way the ferrite itself is specified by 
a single parameter, its saturation magnetization, and its state by an- 
other, namely the applied field. The object of the calculation, then, is 
to find, for a given ferrite and a given guide radius, the mode spectrum 
of the wave guide and the variation of propagation constant with mag- 
netic field. This is done by a semi-graphical method supplemented by 
exact analytic foi'mulae in the neighborhood of certain critical points, 
series expansions in cei'tain regions and some numerical computations 
in others. A sketch of a similar procedure applicable to the plasma is 
given. 

It should be pointed out that the filled cylindrical waveguide is not a 
topic of the highest importance from the technical standpoint. It is for 
this reason that no effort is made here to obtain a comprehensive body 
of exact numerical information about the modes. One wishes, on the 
other hand, to exploit the simplifying features of the problem (as con- 
trasted with the more useful case of a cylinder of ferrite not fiUing the 
guide) so that the discussion may be exhaustive, in the sense that the 
complete mode spectrum is exhibited. 

In Part II we deal with cases of transverse magnetization. By that 
term we mean the following: the microwave fields propagate in a direc- 
tion normal to the dc magnetization and they do not vary along the 
magnetization direction. They may then be separated into two inde- 
pendent sets of field components, of which only one explicitly depends 
on the dc magnetizing field. For these two fields wave impedances are 
defined which can be used for matching purposes. A few simple examples 
are then given. One special case, that of the "non-reciprocal helix" utiliz- 
ing ferrite, is of importance in traveling-wave tube work and is discussed 
at length.' The slow-wave propagation along both a cylindrical and a 
"plane" helix are treated; magnetic less is analyzed in some detail for 
the plane case, and general rules are given for its approximate deter- 
mination in the cylindrical case. 

In Part III perturbation theory and some miscellaneous topics are 
taken up. Suitable perturbation methods are developed for cases in 
which the wave guide fields are drastically modified over small volumes 
(as occurs if thin pencils or thin tliscs are inserted) and also for situations 
in which the local properties of the medium arc but slightly disturbed 
over finite volumes. Among the miscellaneous topics dicusssed is the 
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propagation between infinite parallel planes filled w ith ferrite in a longi- 
tudinal magnetic field. The effect upon Faraday rotation of multiple 
reflections is considered. 

2. THE PHYSICAL PROPERTIES 

The propagation of electromagnetic waves in a medium is governed 
by Maxtt'ell's equations which connect the space variations of E and H, 
the electric and magnetic intensities with the time variations of D and B, 
the electric displacement and magnetic induction. To characterize the 
particular medium relations may be given of the form T) = \\ e\\ E and 
B = II /i II ^ where I| c || and || m || are the dielectric and permeability ten- 
sors. For disturbances whose amplitude is in some appropriate sense 
small, the elements of these tensors mil be independent of rf ampli- 
tude, but mil depend upon the dc state of the medium, upon the fre- 
quency of the signal and in unfavorable cases upon the wavelength of 
the latter. With the assumptions made in this paper the dependence 
upon wavelength mil not arise. 

The form, of || e 1 1 and 1 1 At 1 1 may be known experimentally or it may be 
deduced from some molecular model of the medium. If the equations of 
motion of the parts of the medium are known under applied electric and 
magnetic fields, the displacement and magnetic induction resulting from 
this motion may be found exphcitly. In isotropic media and in the ab- 
sence of applied dc fields, each component of the displacement or of 
induction depends in the same way upon the associated component of 
E or H. The tensors then become diagonal with equal elements. The ap- 
plication of a dc magnetic field, say in the 2-direction, causes ions to circle 
about this field or magnetic dipoles to precess about it. It follows that a 
rf electric field in the ionic case or magnetic field in the ferrite, normal to 
the dc magnetic field, will produce a component of motion at right angles 
to itself and in time quadrature mth it. From symmetry and from the 
equations of motion in a magnetic field the tensors may be expected to 
be now of the form 



(1) 



where a is an even function of magnetic field and h an odd function, c, m 
general, will be mdependent of the magnetic field. 

That a and b at a given frequency and for a given sample of the medium 
are not independent but are related through the magnetizing dc field, Hq , 
is a fact of which we need not take cognizance when solving Maxwell's 
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equations subject to the appropriate boundary conditions. Their solution 
will determine the propagation constant ^ of a wave as a function of a 
and h, no matter what their interrelation. On the other hand, in a given 
experiment /3 is generally' determined as a function of one parameter 
only: the magnetizing field //o . Comparison of the family of calculated 
results ^ = &{a, h), with the results /3 = ^(Ho), found experimentally 
will, of course, determine a and h as functions of Ho . 

If, howe\'er, A\"e have a prior kno^^ledge of a and h in terms of Ho , 
either through postulating the correct dynamical model for the medium, 
or through independent experiments, "sve can utilize the functional form 
of a and b in our analysis of ^, and thus arrive directly at ^ as a function 
of //o . 'f'he distinction between the two methods is by no means aca- 
demic; early introduction of such a functional form of a and b into the 
waveguide problem actually simplifies the analysis. Aside from this prag- 
matic consideration the latter method seems to us more appropriate for 
another reason: it is hardly the task of analysis of technical devices to 
check on the physical theories that give a and h as functions of //o ; such 
checks are made by experiments specifically designed to avoid the ana- 
lytic complexities attending the solutions for most of the technically 
important structures. 

Accordingly we adopt the more direct approach of expressing a and h 
in terms of //o (and, of course, in terms of the magnetic or electric carrier 
density of a given sample) throughout these papers, even in those few 
cases in which ^ can be expressed analytically as a function of a and h. 

2.1 Ferrites 

Most ferrites used in microwave applications are fully saturated in dc 
magnetic fields that are small compared with the dc field with which they 
are biased in operation. We shall therefore always postulate a fully 
saturated sample. Accordingly the magnetization vector M at a, point in 
the sample will always be of constant magnitude, although its orienta- 
tion will change in the ac field. 

One equation of motion for M that takes this into account is 

^ = y[M X Ut] - r^ [M X 11/ X Ht]] (2) 

at \ M \ 

where Ht is a total efi'ective magnetic field seen by the spins that make 
up M, I is the time and 7 is the gyromagnetic ratio appropriate to elec- 
tron spins, whose f/-factor is close to 2. The expression on the right hand 
side of (2) is in the nature of a torciue; the force on M is always at right 
angles to M) thus leaving its magnitude unchanged. The first term on 
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the right of (2) is quite well substantiated by quantum mechanical con- 
siderations. It is a vector normal to M and to the force Ht and is re- 
sponsible for the precession. The second term is also a vector normal to 
M, but is in the plane of M and H in a sense such as to reduce the angle of 
the precession. It thus represents a damping. Not much is known about 
the precise mechanism of the damping, so that its phenomenological 
representation by the second tenn of (2) is still in doubt. 

Ht , the total field acting on the electron spins, is made up of terms 
not all of which are of electromagnetic origin. It consists of the dc field 
Ho mthin the sample, the ac field H, the anisotropy field, and the field 
ascribed to the quantum mechanical exchange forces between spins. 

^0 in the sample must be calculated from the applied dc field ^ext by 
a purely maguetostatic calculation, which, in the case of sufficiently 
simple shapes, can be carried out mth the help of the appropriate de- 
magnetizing factors. Throughout this paper it is assumed that this 
problem has been solved, so that Ho is given. Furthermore it is assumed 
that ^ext and Ho are unifonn. Boundary effects due to non-uniformities 
of Ho are neglected. 

The microwave field H in the sample is one of the unknowns of the 
problem of propagation, and mil appear in the solution of Maxwell's 
equations subject to the appropriate boundary conditions. 

The anisotropy field, a property of a single crystal of ferrite, arises 
from the fact that through the medium of spin-orbit interaction, the 
electron spins can "see" the orbital wave-functions. Since these have the 
symmetry properties of the crystal, it is to be expected that the aniso- 
tropy field will be a vector function of M, with the symmetry properties 
of the crystal. The samples of ferrite used in practice contain a great 
many small crystals randomly oriented, so that the net effect of the 
anisotropy field on microwave propagation must be obtained by means 
of an averaging procedure. The integrations involved are laborious and 
have not been carried out so far. We shall therefore neglect anisotropy 
altogether. Since anisotropy fields are usually of the order of a few 
hundred gauss, this wiU put our results in error below fretiuencles of 
about 3,000 mc/sec. (Corresponding to a precession frequency of yHo = 
3,000 mc/sec, Ho is about 1,100 gauss.) 

The field between two spms ascribable to exchange forces will be zero 
when the two are parallel, and thus arises out of differences of^spin ori- 
entation (that is, differences of M) from place to place. In fact, analysis 
shows that this magnetic field is proportional to V^M for cubic crystals. 
Thus equation (2) really involves position coordinates as well as time. 
Hence the ac part m of M at a point mil depend not only on the ac field 
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H at that point, but on values of H throughout the volume of the sample. 
Therefore B, whifh is nolj + m, will likewise be a functional of H over 
the whole sample. Fortunately it turns out that the spatial variation of 
^ in a microwave structure is so much slower than that characteristic of 
the "spin waves" to which VA/ gives rise that this effect is quite negli- 
gible at microwave frequencies. Only in the most immecliate vicinity of 
gyromagnetic resonance could such effects become significant. 

Thus, we shall regard Ht simply as the sum of the dc and ac magnetic 
fields, ^0 + H, and correspondingly M as the sum of the dc magnetiza- 
tion (directed along Ho in a saturated sample when anisotropy is neg- 
lected) plus an ac part m. Equation (2) must now be solved for in in 
terms of H. It is a non-linear equation, whose solution m will depend on 
H non-Iinearly, as will B. Even if m could be determined in this way, 
Maxwell's equations would become non-linear, and hope of their solu- 
tion remote. It is therefore necessary, and in the great majority of ap- 
plications also quite sufficient, to assume that the ac quantities in (2) 
are so small that their products can be neglected and only linear terms 
taken into account. The terms m and H may now be assumed to vary 
as exp jo}t. 

Under these circumstances, (2) becomes 

^ = 7{[m X Ih] + Wo X //]) 
at 

""^ (iMo X [m X Ho]] + [Mo X [Mo X IJ]]), 



M, 



and is easily solved for m in terms of H, and of the dc quantities Ho , Mo 
which we shall as.sume to point in the z-direction. Each of the components 
nil , my is a linear function of both H^ and Hy and when they are sub- 
stituted in the components of the eciuation B = y.oH + m, lead to ex- 
pressions of the form (1) for B in terms oi H: 

Br = m//i - jnHy , 

By = jkH^ -f mH„ , and (3) 

B: = iioH^ . 

It is convenient to introduce two auxiUary ciuantities 
I 7 I Ho I 7 I Mo 

CO tlodi 
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and in terms of these one obtains the relations first derived by Polder: 



MD ^ (t\1 + «^) - 1 + 2jaa Sgn p 



(4) 



K — V 



where the function 

sgn p = + 1 p > 

- -1 p < 

IT is the ratio of the natural precession frequency ~ ] 7 [ Ho to the 

signal frequency, p is the ratio of a frequency — | 7 1 Mo/iio , associated 

with the saturation magnetization Mq , to the signal frequency. Note 
that a and p always have similar signs; if Ho is reversed, so is the satura- 
tion magnetization. Equations (4) are true only for a fully saturated 
sample. Therefore they hold good only for values of a greater than the 
very small value corresponding to the amount of Ha required to saturate 
the sample. In practice that value of Hq is generally so small that this 
restriction is trivial. In the text a number of formulae will appear 
which apply "near 0- = 0". These are to be understood as applying near 
the very small value of a that corresponds to saturation. 

Equation (4) has an interesting implication with regard to the loss 
parameter a. If a were zero, we would have 

^ -1- .^^„, and 

(5) 



MO 1 - ff^' 

K P 



HO 



1 - <7^' 



and these equations describe the loss-free case. If in equations (5), a is 
replaced by (a -j- ja sgn p), the resulting expressions check (4) to order a. 
For small a, it follows that any propagation problem need be considered 
for the loss-free case (5) only.* The first order change due to loss in any 
formula so obtained can be deduced by differentiation of the formula 

* A form of the damping term in Equation (2) , no less juetified experimentally 
than the one used above, is - r-TT, ( ^ X 'T^ )■ ^^^" *'^'^ expression is used the 
permeabilities are exactly functions of the variable, a -\- ja sgn p. 
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Figs. l(ii) and 1(b) — Tlie rehitive permeabilities n/^o and k/ho versus <t. 



with respect to <t, and multiplication hyja sgn p. Of course this procedure 
is invalid close to resonance (o- = 1), when terms in a play an important 
part. Equations (5), which mil hereafter be cahed the Polder equations, 
are plotted in Fig. 1. The quantity 



H 1 — pa- — ""^ 



G" 



is shown in Fig. 1(c). It occurs in the waveguide theory, and also in the 
theory of other microwave circuits considered later on. The ratio m/mo = 
y-fnn will be denoted by j-n . At a fixed p, fi/iiQ decreases from unity at 
a- = 0, through zero at a ^ —p/2 -\- ■\/pV4 + 1 to ~ co at o- = 1 — 0, 
and then from +coatff = l+0 steadily down to unity at tr = =© . k/^o 
increases from p at o- = to + <» at tr ^ 1 — 0, and then again from 
— ^ at (7^ 1+0 to zero at o- ^ ^ . 

It has already been mentioned that the anisotropy fields are of the 
order of a few hundred gauss. For most ferrites the saturation raagnetiza- 



588 THE BELL SYSTEM TECHNICAL JOUENAL, MAY 1954 

tion is about 1,000 gauss. It -\^dll therefore be consistent with the neglect 
of anisotropy to assume that the appUed frequency is such that p is less 
than unity and this will be done hereafter. 

2.2 Ion clouds or plasmas 

Since these are considered in much less detail in these papers, their 
physical properties are stated only briefly here. 

Instead of a tensor relationship between B and, H we now have one 
between the displacement vector Q and the electric field E. If the mag- 
netizing field is along the z axis, we have 

D, - €E, - j-qEy , 

Dy = jriE^ + eEy , and (6) 

D. - e^E, . 

If the medium consists of equal densities R of positive ions and elec- 
trons, and if collisions and thermal velocities are neglected, e and rj can 
be calculated for weak ac disturbances ^e""' from the equation of motion 



h = _1 Ec''-' + y[v X IJ], 



in 



where v is the velocity vector of the electron and y ^ efio/m, in the usual 
notation. When this equation is solved and the abbreviations 

too = r -Ho; ff = — ; q = — , "p ~ ,^ 

are introduced, one obtains, from the fact that the total current is 
joieoE -\- Rv, the heavy ions being assumed stationary, 

^ = ^0 f 1 + 



cr^ - 1 



.2 



J? = fo 



qa 



, and (7) 



c^ - 1 
e, = eo(l - q^), 

where to is the dielectric constant of vacuum. The waveguide theory will 
involve the parameters 

2 2 

v^ = t = I ^ — ^1- , and 



(1 - a^)(L - q 



(8) 



_ V _ q '^ 

"^ - 7 - a^ + 9^-1 
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These results apply to stationary plasmas only. If the plasma were 
an electron stream moving along the wave-propagation direction, for 
example, the dielectric constants would depend on wavelength also, and 
the propagation problem would be much more involved. 

The variations of e and tj with a are shown in Fig. 2. e/co for a given (/ 
starts at 0- = 0, € = cd {1 — q'), decreases through zero at tr = \/l — q^ 
to — 00 at cr = 1 — 0, starts again from + oo at c = 1 + 0, and de- 
creases to Co at 0- = 00 . 7j = when a = 0, decreases to — ooatff= 1 — 
and then decreases from + » at 1 + to zero at tr = cc . 

We note that similar formulae apply to the electron-gas in semicon- 
ductors at temperatures sufficiently low and frequencies sufficiently high 
so that damping is not important. However, the formulae have to be 
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Fig, 1 (c) — The rutio p^ = ic/ti versus o-. 
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generalized in some of those cases to take into account the existence of 
groups of electrons with different "effective masses" m, 

3. THE SOLUTION OF MAXWELL's EQUATIONS 

Maxwell's equations will now be solved in a cylindrical waveguide 
filled with a hypothetical medium which contains the ferrite and the 
plasma aw a special case. It will be supposed therefore that both its per- 
meability and dielectric constaut are tensors of the form previously 
considered. 

3.1 Field components 

The following notation va\l be found convenient. The projection of a 
vector A upon the plane normal to the z-axis \A'ill be written At . If the 
components of At are a, /3 then an associated vector having components 
(/3, —a) is devoted by A,*. A similar notation is used for differential 
operators. Thus, if V denotes {d/dx, d/dy), V* denotes (d/dy, - d/dx)t. 
Denoting scalar products by a dot, the following identities are evident 

At*-At* = ArA, ; (At*)* = -At ; AfAt* = 0; 
ArBt* = -At*-Bt\ 
and 

ArBt* = 2-component of \A X B]- 

Also if /; is a unit vector along the positive 2-axis, k X A = —At*. 
Similar relations hold for differential operators. If one denotes the star- 
ring operation hy the symbol P then clearly 

p' = -1; P-' = -P; -^ = ^4-2 (« - ^)' 

P -\- a 1 + a' 

where a is a number. 

Max\\'eirs equations may now be written, for that case in which the 
dependence of any component upon t and z is of the form e' "'^^^ , in the 
form: 

V-Ht* = jW,B., 
V*£. + j^Et* = -joiiiHt - ukHi*, and 
V-Et* = -jc^u^M., 
where use is made of equations (3) and (6). 

t The operator V* is called "flux" by Schelkunoff. Strictly, one should write 
V( and V(*, rather than V and V*, but this is needlessly cumbersome. 
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Figs. 2(a) and 2(h) — The relative dielectric constants e/eu iiikI ij/«u versus <r. 



It is desirable to remove scale factors as far as possible. A unit of 
length given by 

1 ^ 1_ 

w'ill be used to measure lengths. This unit is Xo/2t, where Xo is the wave- 
length in an unbounded, unmagnetized medium. It T\ill be assumed that 
/3 is in future measured in units of (3o ■ Finally all magnetic fields will 
be multiplied by x^n^/tr to give them the dimensions of electric fields. 
Vsing the definitions of the p's and p's given in Section 2, Maxwell's equa- 
tions may be put into the form: 






(10a) 
(10b) 
(10c) 



and 



V-^,* = -iT/, . 



(lOd) 



592 THE BELL SYSTEM TECHNICAL JOURNAL, MAY 1954 

Et and //( may now be eliminated yielding two simultaneous second 
order eciuations for E, and H, . These, in turn, may be combined to 
produce two independent second order equations each of which is satisfied 
by an appropriate linear combination of E^ and //^ . These equations 
may be solved and E, and H, expressed as linear combinations of the 
solutions. The transverse fields are then written in terms of E^ and H^ 
and, finally, the boundary conditions are applied leaving a transcendental 
e(iuation in /3^ 

Operating on (10a) and (10c) "withV- and taking account of (10b), 
(lOd), one finds that 

j0S7-Ht* = veU'^-Ei+jpeH,) = -^i^\,and 

3fiV-Ei* = -VH{N-Hi + jphE,) = 0H, 

Operating on (10a) and (10c) withV*-, usingV*-V* = V'' and so on, 
one obtains, using (lOb) and (lOd), 

V'H. + j0V-Hi = PE(-Ih + peSZ-Ei), and 

V%+j^V-Et = -i'n(E,-\- PHS/-Ht). 

Now, elimination of V -^i and V -Hi between (11) and (12) yields 



V'H, + rs ( 1 - pe' - —J H. = Mpe + ph)E, , and 
V'E, + vJl- pi' - —)e. = -Mpe + P>f)H., 



(13) 



equations which demonstrate that pure TE or TM fields no longer exist, 
as the result of the presence of p's. H^ or E^ might now be ehminated 
between these equations giving a single equation in V" and (V")^ but it 
is more convenient to find those linear combinations of E^ and Hz which 
satisfy a first order equation in V^ Writing such a linear combination as 

,/- = £,+ iAff. , (14) 

and adding jK times the first of equations (13) to the second, it is found 
that this is an etiuation in ^p alone of the form 

VV + xV = 0, (15) 

provided that A is a root of the quadratic 

Ve \l - PE - — } - Ph [I - Pn - I 

A^ - _X '-^ ^ "^^^ A - 1 = 0. (16) 

Kpe + Pu) 
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2 



The value of x is then given by 



= - JiU 



Xi,2" = VE [l - PE- - - 0(pe + ptf)A2.i , {17a) 



or 



X1.2' = Vh(i - PfT - — ) + Kpe + Ph)M.2 , (17b) 

where Ai and A2 are the roots of (16) and xi, X2 are the corresponding x'- 
The labelling of the roots is not important, but consistency must be 
maintained, From (14) E^ and H^ must satisfy 

E, + jAiH, = Vi , 
and 

E, + JA2H, = ^2 
so that 



and 



p A2^1 - Ai^2 ,^o , 

-^j = —, r— , (18a) 

A2 — Ai 



H, = j p^ . (18h) 

A2 — Ai 



Solutions of (15) may now be sought in cyhndrical coordinates. To 
satisfy the boundary conditions in circular guide it will be necessary to 
assume the solutions to vary as f"'*, where 6 is the polar angle and n is 
any integer, positive, negative or zero. Equation (15) then becomes 

if r is the radius. Solutions which are regidar mthin the guide will have 
the form of constant multiples of /^(xi.sr), where J„ is the «"' order 
Bessel function. The solutions of (15) are, then, 

^^1.2 = -■li,2J„(xi.2/-)e'"'', (19) 

where the ^'s are constants. Ez and Uz can be found now from (18), but 
further equations must be found to express Et and Ui • Usuig P to denote 
the starring operation, (10a) and (10c) may be re-written as 

U*'fP - peve)Ei + jmt = -VH,, 
iuid 



....; 
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-j^Et + (jvsP - PBVn)Ut= VE,, 

which yield 

{[veVhO- + PePh) — (3^1 — j^HVeipB + pE)P\Et 

= -j^VE, - vhUP - Pi,)^H,, 
and 
{\vbVh{1 + PbPb) — /3'] - JvhVe{ph + PB)P\Ht 

= vbUP - pe)VE, - j^VH, 

The term in parentheses may be removed by using the rule for inverting 
such expressions in P which was given earlier. This process gives 

m = [(l + PEPB - ^^ + J(PB + Pe)p] f20^5 

[- j^VE, - vbUP - pb)VHI 



and 



nUt = 



( 1 + PePh - -^ J + J(PB + Pe)P 



(20b) 



UeUP - pe)VE, - j^VHl 
where 

S2 = VeVh ( 1 + PbPb — ■ } — \Pe -\- Pb) , 

\_\ VeVb/ 

= y^y^ \A. - (1 + Pe) (1 + p^) J- ^ (1 - p«)(l - p„) • 

\_yEVu J L'^Ei'tf J 

It may be noted that for plane waves in the unbounded medium along 
the z axis, which have E^ ^ Ih = 0, £2 must vanish and that the propaga- 
tion constants for such plane waves are evidently given by 

/5" = VeVb{1 ± Pe) (1 ± Ph). (21) 

The values of E, and H^ given by (18) may now be substituted in (20) 
and the operator P removed. This gives, finally, 

(Ai - k2)^Ei = j\ f 1 + PePit —j (1SA2 - PhVb) + viiipi, + Ph) Vl^i 

r ^"/ &' M ^^^^' 

- (/5A2 - PHl'ff)(pE + pff) + I-ff ( 1 + PEpH - I V*!/'! 

|_ \ VeVii/ _ 

minus the same expression mth suffixes 1 and 2 interchanged. 
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(Ai — A2)S2Wi = ( 1 + pspn — ) {A-ii'BpE — (3) — VE^.2{pE + Pn) V^^i 

r / &'\ 1 ^^^''^ 

- j Vei\i I 1 + PePh ) - (pfi + Ph){A-iVhPe — 0) V*^i 

L \ VeVii/ J 

minus the same expression with suffixes 1 and 2 interchanged. 

Equations {22a) and (22b) may be written in a variety of equivalent 
forms by making use of the relations between Ai and A2 . The manipula- 
tions which have been used in deri\'ing (22a) and (22b) assume the use 
of rectangular coordinates, but the results are vaUd in polar coordinates 

*» 1 3 \ 

if El means {Er , Eg) and V means ( z- , - t:;; ) • That this is the case may 

^ or r 06 

be seen from the consideration that the rotation, —6, which carries the 



vector (E^ , E«) into the vector (Er , Eg) also transforms I — , — ) into 

dx dy/ 

d_ Id^ 

dr * r d9^ 

3.2 The characteristic equation 

The boxuidary conditions of the problem are that E, = and Eg = 
at r = ro , the radius nf the guide. E^ is given by [see (18)1 - 

(A, - AO/i', = [A,.4,y,.(xir) - A,A,J„ix2r)W''\ (23) 

and vanishes at r = ro if 

. ^n(X2'*o) . ■/n(Xl''o) 

Ai = 7-^ ; A.2 T . 

A2 Aj 

Hence the relations hold: 

A2.1 



From (22a) it follows that 

(Ai - A2)ilEs - 



./n(x2ro)c^'"' r nf/, , 0' 



1 + PePh — ) (/3A2 ~ pnys) 

r l\ veVh/ 



+ vnipir + Pk) j* JnCxi'') + U/SAa — PhVh){pe + Pii) 

+ Vh\\- -V PePh — 

\ VeVn 



i\fn'{x\r) 



* In A])pen(iix III tlic iicld ooinponents in polar coordinates are written out 
fully for the feirile :iiitl phisina csises with some changes in notation which arc 
introduced in Sections 4,11 and 4.2. 



>x .-.- - 
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miuus the same expression with the suffixes interchanged. Hence 

(Ai - Ao)SiA'fl{ro) 

XVeVh / 

f W XiroJn'(xiro) 



To 



(24) 



+ |/3(.. + ,«) - A... \X - pu - — )| -^^-^^^ 

where use has been made of the relation A1A2 = —1. Therefore the 
characteristic equation for ^ is obtained by equating the term in square 
brackets to zero. Because of the quadratic relation satisfied by A and 
the relation between A and x, it is possible to write the characteristic 
equation in a great variety of ways. It mil be convenient to introduce a 
function 

FM ^ F..{x) = ^^ . (25) 

Using the F-function and replacing the A's by x's the characteristic 
equation may be written :* 

nvsixi — X\) 



_VeVh 



P£(l — PH ) 



1 =^V„(xiro) 



/3(P. + P.) A.^ ^^^^ 

- ^ f„(X2ro). 
Ai 



The asymmetry of this equation between pa , vh and ps , "b arises from 
the fact that the boundary conditions involve electric field components 

alone. 

It may be noted that if the basic solution had been taken to vary as 
cos nd or siti nB, the expression for Ee would have been a linear combina- 
tion of sin nd and cos nB that could not have vanished at the walls for 

alU. 

In passing we remark that for a guide of arbitrary cross-section, the 

* The cliaractcristic equations given in Reference 4 were apeeinlizations to the 
ferrite and |)l!isma. crises of the form in square hnu:kots. They have also been de- 
rived hy Kales^ and Gamo^. These authors have given expressions for some, though 
not all," of the varieties of cut-off point derived in this ])ii|)er and classified them 
as TE or TM according to the field configuration at cut-off. By contrast, they are 
classified here by their association with quasi-TE or quasi-TAI limit modes which 
reduce to the usual TE and TM modes in the uumagnetized medium. 
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boundary \-alue problem may be put into the form, of which (26) is a 
special case, 

and 

V'/2 + X2V2 = 0, 

. /X^ ^ _ Xl df2 \ ^ Vh{ x2' - Xi) \PH^ 4. n - ''\\ ^ 
^ \A, dN A. ON J 0Cp. + P») ^eph '' ''" 'i dS' 

where d/dN and d/dS are normal and tangential derivatives at the guide 
surface, where, in addition, /i = /2 . 

4. DISCUSSION OF THE PROPAGATION CONSTANTS 

At this point we specialize the characteristic equation (26) to one or 
other of the two media. 

4.1. Theferrite (ps ^ 0, ve = 1) 

4.11. After some rearrangement the characteristic equation becomes 

1 ..V..MW ^2^^ 



1 [ F.ixi. 



.(xiro) 

— — ?i 

xi' ■ 



FniXiro) 

— ^ -- n 



where X2,i = /3 Ai.a and the X satisfy 



(1 — J'h) I 1 — — ) + VupH^ 

X:./ - ^ '-^ X,.2 - ^^ = 0. 

PH 

The x's are given by 

2 /, /5' 



(28) 



Xiy = (1 - ^1 -ph\i.2. (29) 

From Polder's equations for pg and ch , (28) may be written 

Xi,/ - [p + ff(l - ^')]X,.2 - ^' = 0, (30) 

or 

XiXa = -P\ (31a) 

Xi + X. = p + (r(l - /3'), 

= p-\- a -\- (rXiXa . (31b) 

If ^ be eliminated between equations (28) and (29), xi,2^ may be ex- 
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pressed solely in terms of Xi,2 , Ph and vh in the form 

" 1 - Xl,2'' 

X1.2' = " \ ~ ■ 

— - 1 

1 - ^^ Xi.. 

PB 

Again using Polder's formulae, this becomes 



2 1 — Xl,2 



x,/ = f^^- (32) 

1 — O'Al.2 



With these expressions for the x, the characteristic equation takes the 
form 



G'(Xi,cr, ro) = G(X2,cr, ro), (33) 



where 



G(X. ", n) = ^^^ 



M-'V^ 



<7X 



— n 



(34) 



Equations (31b) and (33) may now be considered for a fixed a and p 
as determining associated pairs of values for Xi and X2 . Such a pair in 
turn determines /3^ = — XiX2 . Since /3^ must be positive for propagation 
Xi , X2 must have opposite signs. The convention will be adopted that Xi 
is positive and X2 is negative. Equation (31b) will hereafter be called the 
Polder relation and Equation (33) the G-equation. 

An important fact of which frequent use \viU be made is that the 
transformation 

Xi -^ — X2 , X2 — > — Xi , cr -^ — cr, p —> —p 

leaves the Polder relation and /3^ unchanged and converts n to —n in the 
G-equation. It follows that it is necessary to consider positive n only, 
provided we allow the pair <t, p to take on negative as well as positive 
values. This corresponds to the physical fact that a right-circular wave 
in a backward-directed magnetizing field behaves like a left-circular 
wave in a forward field. 

The discussion in this paper is confined to the first azimuthal mode 
number n = ±1. Accordingly the symbol F will replace Fi in what 
follows. 

Before commencing the graphical analysis of the G function it is ad- 
vantageous to consider briefly the function /^(.r) = xJi'(x)/Ji(x), which 
we require for real and for purely imaginary x. By logarithmic differentia- 
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Fig. 3 — The function Fix) = 



J,(x) 



tion of the infinite product for Ji(x), F(x) is found to be given by 



2 ^2 



where the^n's are the zeros of Ji(x). 

Thus, F{x) is real if x'^ is, which is always the case here. For positive x, 
F(x) is an always decreasing function of x, which has an infinite number 
of first order zeros and poles. The zeros are those of Ji'(.r) and will be 
denoted by u,, . The poles are the zeros of Ji{x). It may be recalled from 
the properties of Bessel functions that for large n these zeros and poles 
are essentially equally spaced with a separation ir/2. When a; is a pure 
imaKinary, equal to Jy, F{x') becomes ijI-i{D)/Ix{tj). This is a steadily in- 
creasing function of y, always positive, and behaving like y — ?^ for 
large y. The function F is shown in Fig. 3. Further formulae pertaining 
to F arc given in Appendix I. The inverse function F~ (.r), which is also 
of some importance, is a multivalued function of .r, whose behavior is 
readily understood from the figure for F{x). We are now ready to proceed 
with the graphical analysis of the C-equation. 

In a rectangular coordinate system with X as abscissa and a as ordinate, 
a contour map is sketched of the function 
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G = 



1 - crX 
1 - V- 



\nn 



1 - X^ 
1 - a\ 



- 1 



for all values of X, a from — « to + °c , ro being kept fixed. This can 

be done as accurately as desired by first drawing the contours x = 

1 — X^ 

= constant (Fig. 4), along each of which G simply behaves like 

1 — (TA 

A/\ + B and is easily evaluated mth the help of a table of F. However, 




-1.5 -1.0 



-0.5 



0.5 

\ 

1 - \' 

Fig, 4 — Curves of constant x^ = i° t"6 a — \ plane. 

1 — ff\ 
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many features of the (?-contours are already detennined by the position 
of the contours Cr = ^ , and G = 0, across which G changes sign (from 
± X to ^ » or from ±0 to =F0). Because of their special role in the 
subsequent analysis it is desirable to introduce a scheme for their enum- 
eration. The infinity and zero curves in the right-hand half-plane will he 
denoted by / and 0, respectively, those in the left-hand half-plane by /' 
and 0'. All but two of the /-curves arise from the poles jn of F. Their 
equations are 

:; r = j«"/''o" n = 1, 2, • • " 

V — ah 

Each of these curves has two branches, one in the half-plane X > 0, one 
in X < aiid these are called /„ , /„' respectively. All /„ curves pass 
through X = 1, ff = 1, all //curves pass through X = — 1, c = — l.The 
lines X = 0, X = — 1 are also infinity curves to be denoted by I a , la 
respectively (As X ~^ +1, G tends to a finite value). 
Zero curves of G are given by 



or in a more readily computable form by 

1 ro^d - X^) 



a = — 



X X[F-'(X)? 



(35) 



The branches of F \\) may be labelled according to the scheme: "0" 
for - oo < [F-\\)f < ii'; "1" for if < [/^\X)l' < J2 and so on. The 
?ith branch of F~\X) gives rise to an 0„ curve for X > and to an 0„' 
curve for negative X. All 0,/ curves pass through X = — 1 , cr ^ — 1 ; all 
save one of the 0„ curves pass through X = 1, ^ = 1. The exceptional one, 
seen to be Ou , is associated with the "0" branch of F^'(X) on which 
F"'(l) = 0. For fixed a, G tends to zero as X ^ ^, hence the vertical 
lines X = ± « are also zero curves, to be denoted by 0„ and 0„' respec- 
tively. 

In a sense the two branches of trX = 1 are also zero curves, to be called 
Oc and 0/. 0^ and 0/ are zero curves only when viewed from "one side." 
In the right half-plane, for X<lasffX^-l — and for X > 1 as 
o-X — > 1 + 0, the argument of F tends to infinity and remains real. 
Therefore G passes through all values an indefinite number of times and 
(tX = 1 is a limit line of all contours, G ^ constant. For X < 1 as 
ffX — > 1 -|- and for X > 1 as (tX ^ 1 — 0, the argument of F is 
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Fig. 6 — Qualitative behavior of G{\, a) at large distances from the origin as a 
function of arc tan a/\. r^ is about 2. 



imaginaiy and F tends to infinity. However F 



1 - X^ 



tends to zero 



1 - o-X 
SO that G tends to zero. 

To complete the picture of the G-function given by the form and posi- 
tion of the and / curves it is necessary to sec how it behaves at large 
distances from the origin. This is indicated in Fig. 5 and also by Fig. 6. 
The latter nhows the value of G at large distances as a function of direc- 
tion. In general, along the line ct = cX + rf (c finite), G will tend to — c 
for all d. For c = tq/J^ (which is the slope of the asymptotes to the 7„ 
curves), G again tends to a constant. Now, however, the constant depends 
upon d and assumes all values from — ^ to -1- ^ as a function of d. In 
the first quadrant the sign of ^-ariation of the limiting value of G with 
direction c is opposite to that of its variation mth d near c = To fja ■ 
Consequently local maxima and minima arise as a function of direction 
between successive /„-curves. This suggests the existence of saddle points, 
which may be \-erified directly. In the third quadrant, the dependence of 
G upon c and d does not give such maxima and minima, and indeed no 
saddle points are found there. Finally it is necessary to consider the be- 
havior of G as 0- tends to infinity, while X remains finite, corresponding to 
(1/c) — » 0. If X remains fixed, then for X> — 1,G— >Tx aso--^=too; 
and for X< — 1,G— >±^asff^=tcc.AsX^O, the curves of constant 

G are asymptotic to Xct = ( 1 — -^ ) — -^X, where B goes from — oo to 

+ X with G. Interleaved \\ith these families of curves are the curves 
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G = ± 00 , which are Xtr= (1 ~t^) + 0(^^)- More detailed information 

on these matters mil be found in Appendix II. 

From the G-diagram it would be possible to determine pairs of X-values 
with opposite signs, which, for a definite cr-value satisfy the characteristic 
equation, but, for a given p such pairs would not necessarily satisfy the 
Polder relation (31b). It is necessary to have a procedure which takes 
account of the latter systematically. Such a method may be based upon 
the fact that if, for a and p positive, the Polder relation is solved for Xi 
in terms of X2 it can be thought of as a rather simple mapping of the 
whole X2-quadrant upon a part of the Xi-quadrant (Xi > 0). Similarly for 
(T and p negative there is an analogous mapping of the Xi-quadrant onto 
the X2-quadrant. 

Considering first the case c, p > 0, the Polder relation may be written 
in the forms 

X. = ' + P-^' = i + 1 + iL^lA = T(W. (36) 

1 — (rX2 <T 1 — (TM 

From (36) it may be seen that the curves X2 = const, transform into a 
bundle of hyperbolae passing through the intersection of o- = 1/Xi and 
a = \i — p; that is, through Xio , o"o , where 



tro = -p/2 + j/|'+l , ^10 - V/2 + j/|+l . ■ 

These hyperbolae have the vertical asymptotes Xi = — I/X2 , and intersect 
ff = at Xi ^ p — X2 . For a fixed positive a less than aa , Xi decreases 
from 1/tr to 0- + p as X2 increases from — » to 0, but when a is greater 
than ao , Xi increases from l/a to o- + p imder the same circumstances. 
Thus the whole X2-quadrant is transformed upon that part of the Xi- 
quadrant which hes between the hyperbola Xi = l/tr and the straight 
line Xi = o- + p. It follows that points in the Xi-quadrant which are, for 
a given p, excluded from this region, cannot be the site of acceptable so- 
lutions of the G-equation. 

Since as has already been stated, the Polder relation is unchanged 
by the substitution Xi — > — X2 , X2 — > — Xi , a — > —a, and p to —p, it 
follows that for a and p negative a similar mapping of the Xi-quadrant 
upon part of the Xj-quadrant takes place. The transforms of the lines 
Xi = const, and so forth may easily be found by using these substitutions 
in the foi-mulae already gi\'en. 

Reference to Fig. 1(a) and (b) mil show that dzo-o are the values of a at 
which fi reverses sign. Therefore we may expect ao to play a special role 



GUIDED-WAVE PROPAGATION THROUGH GYROMAGNETIC MEDIA 605 

in the propagation theory, as also does a = 1. The following scheme 
exists: for < o- < ffo , k and fi are both positive; for (tu < o- < 1, k < 
and ti < 0, for a > I, K is negative and ^ positive. If <r is changed to 
— (r,fj. goes into m, ^"d '^' i»to "''''■ 

The procedure which will now be used to discuss the solution of the 
characteristic equation, observing the Polder relations, begins by writing 
the eciuation, for tr, p positive, in the form 

We are already in possession of a contour map of the left hand side 
of this equation in the quadrant cr > 0, X > 0, and of the function 
GQ^i , a, To) in the quadrant A < 0, (j > 0. The latter surface has now to 
be transformed into one in the Xi-quadrant by the relation 

\, = T{\,) = (a + p - X,)/(l - ffX,) 

(or equally well, Xi = r(X2). This may be effected by considering the 
transformation of curves G(X2 , ff, To) = constant, onto the Xi-quadrant. 
For the /' curves whose analytical expression in terms of a- and \- is very 
simple, the corresponding expUcit expression of the transformed curve 
in Xi and ff is simple. Contours other than F are most easily transformed 
by replotting G(Xa , ff, ^o) = const, in the hyperbola-mesh formed by the 
lines T(\2). However, information about particular points and about 
asymptotic behavior of these transformed curves is available in analytic 
form and is stated in Appendix 11. The two surfaces so obtained will in- 
tersect in various curves, along whose projections on the X — p plane both 
Polder relation and G-equation are satisfied. For each such projection Xi 
is a function of a, Xa is then known in terms of a and p, and finally 0' = 
— X1X2 is known. In most cases the general course of the.se curves can be 
found mthout resort to much numerical analysis. Each of the curves is 
associated with a definite mode and it follows that the classification of 
the modes can be carried out fairly easily. The approximate location of 
the solution curves relies upon the fact that if the position of the infinity 
cui-\'es of both surfaces is kno\\Ti, continuity considerations will fre- 
quently assure the existence of an intersection within certain regions. 
Moreover, the neighborhood of (certain special points on these solution 
curves can be investigated analytically. These are points at which one 
or both of the G-fimctions may be approximated by a simpler expression; 
included among these is the point at infinity. 

It is clear that for a and p negative the whole procedure outlined above 
may be carried out in a similar way, with the tr > 0, X > quadrant now 
being transformed on to the ff < 0, X < Cfuadrant. 
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It is possible to translate such solution curves into jS^ — a curves in a 
direct graphical manner if a mesh of constant /3^ lines is drawn in the 
first quadrant. From (30) these are given by 



K'- [p + <7{l - 0')]X,- 0' = 0, 



or 



1 " Xi^ 



1 - 



a + 



-^ = 1 - ^^ 
Xi 



1 - ,32J 

The contour, = b,is just the contour x^ = I — b"^ displaced along the 
o--axis by an amount, — p/(l — &^). The contours of constant /3 all pass 
through the point a = ao , X ^ ao -\- p and are shown in Fig. 7. Their 
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Fig, 7 — Contours of constant /3^ = X for p = 0.5, 

1 — ffX 

Ij9»(X,a, ;») =^=(-X, -a, -p)j 
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course in the third quadrant is immediately found by refleetion in the 
origin. 

When p = the magnetization of the ferrite vanishes and it is clear 
that we should then obtain just the modes of a guide filled with isotropic 
material (p. = ti^ ; k ^ 0). Superficially it might appear that, sinc^ the 
equations {31b) and (33) depend upon a, even for p = 0, this result 
might not be attained. We now show that ^^ is indeed independent of a 
for p = 0. It may first be noted that in tfiis case if a s^ 1, the Polder re- 
lation (31b) transforms 



1 - X 



1 - <t\, 
The G-equation reads 
1 — X2a" 



]- into 



1 - Xo 
1 - ffX^ 



and Xi into 



g- - Xo 
1 - <r\z 



£Kv^s-: 



_ 1 - \ia 
" 1 - \.2 

Since Xi 5^ X2 we must have 



1 - aXo 

"o- - X. 



Fir, 



Fir, 



1 - X.^ 
1 - crXn 



= or CO, 



1 - Xr 



' 1- 1 

1 — a\2 



or 



1 - \l.2' 



1 — (rXi,2 



To' 



Thus Xi,2 are roots of 



or else of 



In the first case 






or ^-, . 






-1=0. 



x^-^-'-^x + f^- - 11=0, 



ro- 



To^ 



-X,X2 = ^^= 1 -^, 



and in the second 



ro' 



3n 



_X,X2 = ^^= 1 -■'-,. 



To' 



u^ . ■ I ■ 
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Thus, when a ^ 1, the /3 values, for p = 0, are evidently independent 
of (T and are just those of an isotropic medium. When ff — 1(^ = « = oa), 
p = 0, |3^ is indeterminate and for p small, there is a small region near 
0- = 1, of width '^p, in which /3^ differs appreciably from the isotropic 
value. The convergence of an expansion of /3^ in powers of p[{61) and (62)] 
shows a marked dependence on cr. 

4.12. The scheme of analysis described above will now be illustrated 
in detail by a discussion for a radius ro between Wi and ji , which, if the 
ferrite were unmagnetized, would propagate the 7'^^n-mode alone. 

Figs. 8(c) and 8(d) show the division of the \ — <j plane into regions of 
positive and negative G{\ a) by the various / and curves. A few con- 
tours of constant G are plotted to indicate the behavior of the function 
in more detail. That part of the X — o- plane which is excluded by the 
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2.2 



j'ig^ 9 _ Genmctrical exploration of the solution curves. The permitted areas 
of the X - (T phiiie are divided by the 0, 7, {0)t , {I)t curves into reRions in which 
G(\,<r) and G(Ti\),o) have like or unlike signs. Shaded regions are those of unlike 
sign's. Solution curves (shown schematically hy dotted lines) must lie in regions 
of like signs. Only the first few and I curves are shown. Fig. 9(a) and (b), ro -^ 
H.O; Fig. !)(<■) and' (d),/'o-' 5.0; Fig. 9(e) and (f), ?<= < ;-n < j; ; Fig. 9(g), ro < n, . 
I ;) I < 1, lliroughoul. The horizontal dashed line marks \ <r\ = \ <ro\ . 



Polder relation is indicated, for p = M, by shading. For other p-values 
the straight portion of the boundary of the excluded region is simply 
translated along the X-axis. 

In Fig. 9(a) the allowed region of the first quadrant is shown again, 
together with the transforms (/./)7-, (Oi')r , and (Th')t of the only 
critical curves I/, 0/ and Ib occurring in the second quadrant for the 
present radius. Regions in which G(Xi , ff) and (?(r(Xi), <t) have opposite 
sign arc shaded ; the common signs in the remaining parts of the quadrant 
are as indicated. (In this diagram p is taken to be ^g). 

From the disposition of the surfaces G{\i , a) in the region between d 



-2.2 -2.0 




1.8 2.0 2,2 



Pig. 9Cb) and (c) — See Fig. 9. 
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and /i and of the surface G{ T(\i) , <t) between {OJ)t and (Ib') t it is evident 
that along any contour such as G(Xi , a) - K, G{T{\i), a) mil take on all 
positive values from to « and, in particular, K. Since this is true for 
any K, it follows that the region between Oi , 7i , (Ib)t contains a solu- 
tion curve. Two points on this curve are immediately obvious: the inter- 
sections of (Ib')t , h and the pomt (1, 1) on (0„')t , 0, . The first is the 
intersection of the curves 



X = 1 + 



1 +c' 






1 ^W' 

1 - a\ 



At the point (1, 1), Xa is - ^^ , Xi is unity and /3^ is therefore infinite. 
Armed -with this knowledge we now investigate analytically the behavior 
of /3" near o- = directly from the original G-equation and Polder rela- 
tions. Writing Xi - 1 + ce and ff = 1 + e, (1 - x')/(l - (tX) is to zero 



-2.2 -2.0 -I.S 




- -0.2 



- -o.e 



-0.8 



- -1.2 



- -1.4 



- -1.6 



Fig. 9(d) — See Fig. 9. 
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-22 -20 -1.8 -1.6 -14 -12 -10 -06 -OP -O.-l -OS 




Fig. 9(e) and (f) — See Fig. 9. 
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order just 2c/l + c. Thus, G{Xi , <t) to zero order is 



1 + c 



2c 



2c 



1 +c 



- 1 



But we have, in this case, G(Xs) — 0, so that 

2c zi 



1 + C 7-0^ ' 

where Zi is the smallest non-zero root of F{z) = 1. From the Polder 
relation, the leading term of X2 is —p/(l + c)e, and consequently the 
leading term of )3' is 

2 



P 



v{^-^. 



2u 



(1 + C)€ 



a - 1 



This analysis is readily extended to the next order term, which is stated 
in Section (4.17). 

From analogous considerations concerning the variation of one G- 



-2.4 -2.2 -2.0 -t.a -1.6 -1.4 -1.2 -1.0 -0.8 -0.6 -0.4 -0.2 




Fig. 9(g)— See Fig. 9. 



616 THE BELL SYSTEM TECHNICAL JOURNAL, MAY 1954 

function through all possible values of the other in the region bounded by 
h , (/flOi' ) (0/)t' we deduce that the solution curve just discussed con- 
tinues into that region and persists as u — > co . For, the asymptote of 

(Oi')r is (7 = ~ X, and between it and X = 1, which is the asymptote of 

111 

{Ib')t , G(2'{^i)> '^) takes on all values between and — co ; in particular, 
the hmited range of values assumed by (?(Xi , <t) in this region. The be- 
havior of the solution curve for large a may be deduced by using the 
asymptotic formulae for curves G^(Xi , c) = g and ^{^(Xi) a) = g which 
are given in the appendix. These are 



-g\i - gF 



and 









1 

— , Xi — 


P+2 


\ Tip- 



1 - Wl^_ 

It is clear that j? at a point of intersection is given by — ro fui plus terms 
of order 1/Xi ; substituting this value in the second equation gives the 
solution curve correctly to order 1/Xi in the form. 

2 
^0 -. 

When the solution curve has such a linear as3Tnptote it is convenient 
to calculate )3^ from the formula 

/3^=1+- h terms of order higher than l/tr 

ff a 

which is readily obtained from (30) . In the present case 

^^ = (l^ ^1l\ h + ^^ + higher terms in 1/tr. (38) 

As 0- — > «= , |S tfinds to the value appropriate to the TEn-mode in an. 
isotropic medium (/i — > ^z = /xo , k — »■ as o- — > <» ). Thereby the whole 
solution curve is classified as specifying part of a TEn-limit mode. 

The remaining section of the TEn-limit mode in the upper half -plane 
is again found in the region between (Oi')t and {Ib')t for <7 < cq . Any 
line a = constant < uo cuts these two curves at two values of Xi . As Xi 
varies between these values, ^(^(Xi), a) varies from to — co ; it is, thus, 
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Fig. 10(a) —3- veraus ]) for smtill values of a — the TEu-Hmit mode. 

clearly equal to the finite (negative) G(Xi , tr) somewhere between. This 
situation persi.sts up to a- = ffo — and a sokition curve therefore exists 
between ff = and tr = o-q . It meets (f = for Xi satisfying 



1 - \:- 






i - X. 



-\\ F{r,VT^rx:^ - i], 

2" LA2 J 



and 



Xi + X, = p. 



(39) 



These equations have been solved numerically; the corresponding /3 
= — X1X2 is shown in Fig. 10(a). For ro between Wi and ji a value derived 
for 0^ from the first three terms of an expansion of 0'^ in powers of p, 
equation (61), turns out to be in very good agreement with the numerical 
calculation up to p = 1 , for a = and presumably is good for small ff. 
At ff(, (the point at which m becomes negative), the solution curve is 
"cutott". However, Ihc corresponding ^' is not zero. As cto is approached 
from below (7(Xi , ff ) -^ and so G{\^ , a) tends to zero. Thus, X2 tends to 



618 THE BELL SYSTEM TECHNICAL JOURNAL, MAY 1954 

0.6 



0.5 



0.4 



0.2 









/ 


^— 





To =5,75 




_^ 


^ 






/ 


/ 


















/ 


/ 




^*V 


Si^.OO 










/ 


/ 










\ 






/ 


/ 

1 


/ 


/ 


^^ 


^^ 


»4.2S 


> 


\ 




/Ci 


xi 


'4.825 
/ 


^ 






\ 


V 


/ 





0.2 



0.4 



0.6 



0.8 



-1.0 -0.8^ -0.6 -0.4 -02 

Vo=5.33 p 

Fig. 10(b) ~- j3* versus p for small values of o- — the TMu-limit mode, 
the negative root Xao (unique for the present radius) of 

y 1 - <7oX2 



Fin 



= Xa. 



The associated /3^ is — XsoXio = — — and is shown in Fig. 11(a). The way 

in which /3^ approaches this value as a- — > o-q can be found and is one of 
the more subtle examples of behavior of a mode near a special point. 
Writing a = oo — 5(r, Xi = Xio — 5Xi , we observe that, since ao + p 

— - = 0, the Polder relation in the form 

Co 

1 <T -\- p ~ 1/a 

Ai = - -h - — 

a 1 — 1TX2 

fully determines -^ ; any variation due to 8X2 vanishes at a- = co . 5X2 can 

be determined from the (j-equation. Near crX = 1 — (X > l),G(Xi,o-)ia 
given by 



■0 , / l - o-Xi 
Ci r X:^ - 1 ' 

which near ao , Xio may be written 



Xi 
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-V^a '"" 




da 

Aio -77 — r To 
a\i 

Xio^ - 1 



dG 



The perturbed G(X2 ,a) which (since G(X2o , tq) = 0) is -— 5X2 + OC^tr) 

0A20 

equals the preceeding expression and gives 



/T- ''0 



5X2= -V5^^ 
Am 



Alo 7T — r Co 

OAi 



1 



dXnJ 



Accordingly, 5/3^ = — X^ 6X2 + 0(5(t), a result which shows that /3^ 
tends to its terminal value along the vertical. It is clear analytically and 
graphically that this mode persists as p — > 0, and must be identified with 
the only isotropic mode for this radius, namely TEn . No other branches 
exist below a = a^ , since ^(Xi , a) and G(T(Xi), a) have opposite signs 
except in the region just considered. 

The two solution curves considered so far are not the only ones; in 
fact the infinity of sheets of the surface G(Xi , u) in the region bounded 
by /i , Oc and (lA')r , Fig. 9(a), intersect the transformed sheets G{T(\i)(t) 
in infinitely many more curves. In the blank areas of that region the 
G-funetions have equal sign, and all these areas must be carriers of solu- 
tion curves, since in every one of them every single contour G(Xi , a) = g 
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Fig. 10(c) — (3^ versus p for small values of <r — the TEu-limit mode 
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Fig, 11 — ^^ at the type 2' cutoff as a function of r^ for various p. 2o' cutoff, 
TEii-limit mode; (b), 2\ cutoff, TMu-iiinit mode. The presence of a curve for 
p = is clarified in the text. 
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crosses all contours of G{T{\i), a), in particular G(T{\i), a) = g. AM 
the iidditional solution curves arising in this way start at o- = 1, Xi = 1; 
the rt^"" of them threads its waj' from one blank region to another, first 
through the intei-section of /„+! mth (Ib')t , then through the intersec- 
tion of 0„ with {Oi')r , and finally comes to an end at the intersection of 
I„ with (Ia')t . At the end point (a ^ 1, Xi = 1), X-. and, therefore, j3^ are 
infinite, (just as for the TEu solution curve). At the end point (/„ , 
Ua')t), Xo , and, therefore, 0' are zero. The a and Xi values corresponding 
to the latter are obtained from the equations 



^ 2 -2 

Xln Jn 



1 — tr,|Xin Vq- 



Xi. = a,+ p. (40) 



It is possible to derive the slope d^'/Sa of the (3^ — o- curves at these cut- 
off points. Near cut-off, the infinity h of G(\i , c) is matched by the 
infinity 1/ of G(M , <t). The (^-equation therefore degenerates to 

Writing a — a^ — -rX^ , Xi ^ Xi^ — 77X2 , expansion of the right hand 
side of this equation to order I/X2 furnishes one relation between x and y; 
the Polder equation furnishes another. The two can be solved for x, and 
so, since to first order 

2 ,, ,*'■ — CTn _ .. StT 

88 = — Xi„Xa = Xln = Al — 

XX 

S0^/8(r may be found. It is found that for convenience in computation, 
the results of this calculation are best presented parametrically. Equa- 
tions (4G-8) represent equations (40) and 50'/8c in this way. Fig. 12 
(a) and (b) show the result of some computations. Near a = 1, 0' = 
OS, these added solution curves behave rather like the TEu curve. The 

1 
leading term in the expansion of /3' in powers of _ " is now 

ru- 



ff - 1 

th 



for the solution curve ending at /„ — {Ia')t . Here.2„+i is the (ji + 1) 
root of F(z) = 1, not counting 0. 

It will turn out later that the infinity of solution curves just discussed 
represents an incipient fonn of the whole mode spectmm; the reservoir 
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from which higher modes are drawn as the guide radius is increased. 
That the propagation of modes which for larger guide radii correspond 
to higher TM and TE modes is possible for limited ranges of a might 
be ascribed to the larger ^-values in those ranges, which cause the wave 
to see an effectively larger guide. This explanation is convincing only 
when 0- > 1 . When <to < a < I , m is negative, and the propagation must 
then be the result of an interplay between m and k. In passing we remark 
that we are here dealing with the propagation analog of so-called "shape 
resonances," which physicists sometimes encounter in resonance experi- 
ments on small spheres of ferrite in cavities. 

We now turn to a discussion of the solution curves for tr < which 
lie in the third quadrant. Fig. 9(b) shows the partition of the region 
allowed by the Polder relation (again for v = %) iiit-o positive and 
negative regions by the various /', 0', {I)t and (0)t curves. Regions in 
which G(X2 , a) and G(T(\n), a) have opposite signs are shaded. For 
tr < — (TO , the question whether a given region of like signs is the site 
of a solution curve may, mth one exception, be answered by the same 
type of geometrical argument as used for tr > 0. The singular area is 
that part of the region bounded by Ib and 0^ m which the G-functions 
are both positive. Here both (?(X2 , a) and G{T{\2), a) are zero on 0^ ; 
G(\2 , ff) goes to « on Ib', whereas G(T(\2), tr) is finite throughout the 
region. No intersection can be predicted, then, by the earlier argument. 
It can indeed be shown (for all ro) that there is no such intersection. For, 
in the case y» ^ 0, the solution curves are /„' or 0„' curves as demon- 
strated in Section (4.11). The region under consideration contains no 
such curves, and hence no solution curves. Thus, for p = 0, since G{\2 , o-) 
goes to infinity on Ib', the surface G{T{\2), o") must lie entirely below 
the surface G(\<i , a). Con.sider now^ for fixed a and increasing p, a point 
on the G(T{\^), a) surface whose height remains unchanged. For such a 
point T{\2 , p, 0-) remains fixed and from the Polder relation this means 
an increasingly negative Xa . Since it can be shown that G(X2 , cr) goes 
monotonically from to «> as Xa becomes more negative, it follows that 
G(T(X2), 0-) continues to lie below G(\2 , a) for all p. 

All other regions of common .sign do carry solution curves. That 
corresponding to the TEn-limit mode begins at o- = — 1, Xs = — 1, 
pas,ses through the intersection of (Oo)t and Oo' and per-sists for indefi- 
nitely large it. The asymptotic formula (38), for /3^ at large a also holds 
as (7 — > — OT , if the signs of both a and p are taken to be negative. The 
behavior of /3' near tr ^ — 1 may be found by the same means used at 

a = I. The resulting expression* is to order — :— - essentially the same 

a -]- I, 



* See Section 4.17 for a more exact formula. 
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as the earlier (37) except that the smallest root of F(z) = — 1 replaces 
that of F(z) = 1. The remammg solution curves confined to the region 
bounded by Oo', 0^ , (Ia)t portray the incipient modes already encoun- 
tered in the first quadrant. Their behavior near a = —I also follows 
(58), associated with the higher roots of F(z) = —1. Their end points, 
the intersections of /„' with (Ia)t, are still given by the parametric 
representation (46-8), due regard bemg paid to the signs of o-„ and p. 

The remaining branch of the TEn-limit mode, lying above a = —ao, 
is found in the triangle between {I a) t , u ^ and Ib . Its end points are 
given by (39) with p negative and by the intersection of Ib' with {I„)t 
which is (T = — (1 + p), A2 = — 1, Xi = 0. Thus the cut-off in contrast 
to the analogous branch for o- > 0, is given by /3^ = 0, it = — (1 + p). 
(When p < —1, the branch does not exist at all.) We note that a left- 
circular plane wave is cut off at exactly the same value of a as the TE- 
mode is in this particular case (see, however, the following sections). 
The slope at cut-off is detennined by expanding the G functions near 
their infinities at Z„ and Ib' and utilizing the Polder relation. The slope 
is found to be 

da p(l - F{n)) ■ ^ ■' 

A further solution curve hes in the region between Oc and (Oo)r for 
ff > — ao . It has no analogue in a guide with isotropic material and 
will be discussed later. 

In the discussion of the mode spectrum for radii between Wi and ji 
three distinct types of cut-off point have already been encountered. When 
larger radii are treated it is found that no other types arise.* In Section 
4.17 formulas relevant to the three types are given. An examination 
of the field components in the neighborhood of the cut-off points is of 
some interest. Cut-off points of type one (intersections of /„ and (Ia')t 
or /„' and (Ia)t), at which p^ = 0, have E^ = and the field is of a 
pure TE-type. The medium behaves transversely as though it had a 
permeability, ju — k/h. Although the field is purely TE at cut-off the 
mode terminating at such a point may in the limit of vanishing magne- 
tization be either a TE- or a TM-mode. This impartiality extendsto cut- 
off points of the other types. Cut-off points of type II [(On')r ~ 0^ or 
(On)T — 0„'] occur at 0" = iffo , where fi = and herc/3' does not van- 
ish. In such cases one of the x's is finite and the corresponding contri- 
butions to the field pattern ciuite normal. The other, however, tends 



* Tliere is an exception to this stiitenient. This is the type designated in Sec- 
tion 4.17 as 2oM which cuts off an isolated mode having nn TE or TM analogue. 
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to an infinite imaginary ^'alue and the associated fields are confined 
very closely to the guide walls. The wall currents are very large and 
essentially longitudinal. Type III cut-off points [Is — (Ia)t] at which 
fi = — Khavc/3' = 0, but the fields are not of a purely TE- or TM- 
type. They ^consist essentially of a rotating, transverse, H, which is 
uniform over the guide. The t-omponents H, , Ee and Er are smaller 
by one order of tr — a-cit-ofr and E. , two orders smaller. 

It should he stressed again that, in general, the modes are never of 
pure TE or TM type. Nevertheless, for the sake of brevity, we shall 
refer to them as such; calling them TE-modes or TM-modes according 
to their limit as the magnetization is removed. 

4.13. We now consider the behavior of the modes as a function of 
radius. The reader will be aided by Figs. 8(a) to (e) and 9(a) to (g). In 
preparation for this it is necessary to examme the movement of the 
/„ , /„' and 0„ , On' curves when ro is varied. It will be recalled that the 
equation for the /„ curves and their reflections, //, in the origin, is 

1 A 3 -2 

1 — A ^ ,7n_ 

1 - ffX ~ W-' 

The contours x(^., o") = (^, where x^ = (1 — ?i^)/(l — c^^) have already 
been plotted in Fig. 4. The /„ , /„' curves are among these, and, clearly, 
for a fixed n, the associated x" tlecrea.ses as vo increases. The course of a 
given pair (/„ , /„') may then be seen directly from Fig. 4. The qualita- 
tive behavior of the pair changes radically only when ro passes through 
the value j'„ . Before it does so, 7„ lies, for X between and 1, above tr = \ 
and tends to o- = ^ as X tends to zero. At ro = j„ , the 7„ and /„' curves 
merge into the lines a = \ and X ^ 0. Beyond j„ , /„ lies below u = X 
for X between and 1 and goes to — ^ as X approaches zero. The /„'- 
curve remains, throughout the reflection of /„ in the origin. As ro — > -» , 
I„ tends to the line X = 1, /„' to X = —1. No I„ curves ever enter the 
region X > 1, tr < 0; no /,/ curves enter X < — 1, o- > 0. It is also im- 
portant to relate the In , I ,! curves to the boundaries of the Polder 
regions. /„ curves cut the Polder boundary a = X — p, of the first quad- 
rant in at most one point. As h, increa.ses from to J„ , this point moves 
from 0- = (To to <r = =o . Thereafter, no intersection occurs at fixed -p until 
ro equals jn/'s/i ~ P'\ it here reappears at o- ^ and moves steadily to 
(J- = 1 — p as Td increases indefinitely. The only intersection with the 
other Polder boundary o- = 1/X, is at X = 1, (t = 1, regardless of ro . 
The 0„ , 0„' curves are given by 



- 1 fl 2 1 - XM 
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if the n'-^ branch of F~^{\) is used. Thus, as ro increases, the successive 
curves either all pass through a fixed point (which can only be X = ±1, 
0- = ±1, n > 0) or move steadily up or down without further intersec- 
tion. An On curve starts from a\ = 1 atro = and falls for X < 1, rises 
for X > 1, as J-Q increases. For large X, since <t •-^ jn'/ro^ (X •- 2), n > 0, 
the 0„ and /„ curves move together with a constant separation. Oo is 
singular, since it does uot pass through X, cr = 1 and falls steadily for all 
X; it tends to a = for large X. TheO„' curves rise from trX = 1 at ru = 
for —1 < X, fall for X < —1. They run parallel to h+i for —X very 
large. For small X there is an expansion 

, = (i-:iy j'i + o(x) 

holding for 0„ and for 0„'. This indicates that for ro < w„ , 0„ goes to 
+ oo and On' to — «= for small X, but at ro = «„+! , 0„ and 0/ merge 
momentarily at 

X = 0, 



2ro' 



< 0. 



For larger ro , 0„ goes to — oo and 0„' to + ^o . Since the union of 0„ and 
0„' takes place at a negative cr, it is clear that 0„ curves, unlike 7„ curves, 
may cross the hne cr = twice. Intersections of the 0„ , 0/ curves with 
the Polder boundary are difficult to examine explicitly and this may lead 
to some obscure situations for < | X | < 1. However, for a > tro , 
since On and /„ have a fixed separation for large 1 X |, this pair escape 
intersection with the boundary at the same value of ro , namely jn . 
Similarly On' and /„+i escape together at ro = j^-^x for o- < — (To . 

We shall now examine the effect of varying ro upon the sequence of 
modes when c > o-q . When ro is less than u\ , a case in which the iso- 
tropic medium would not propagate, no part of Oo' hes in the upper 
half plane and there is then no (Oo')r curve. The solution curve which in 
the previous discussion of Section 4.12 was assigned to TEn , after 
passing the intersection [7i — {Ib)t] can no longer escape to infinity and 
terminates on [7i — (Ia)t]- Thus, the TEn mode at this radius has 
become an incipient mode with cut-off and other properties given by the 
formulae already quoted for such modes. As n approaches ?(i from 
below, the ^^ — a curve is double valued between o-cut-ort and some larger 
value. This is borne out by the fact that d^jda becomes positive at 
cut-off, and by the observation that the solution curve bulges towards 
large a between /i — {Jb^t and its terminus. The part of the ^ — u 
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curve along which -j— < 0, ^dll tend smoothly towards the /3' — a curve 

da- 

for To just greater than i(i . The course of the TEu solution curve remains 
qualitatively unchanged for all tq > Hi . 

When ro passes through iti , and the TEn solution curves escapes dis- 
continuously to infinity, the solution ounces below it disengage from 
their former end points 7„+i — {/a)^ and instead end at the point 
I„ — (Ia)t • When ro exceeds ji , tlie curves /i and Oi escape intersection 
with 0" = X — p simultaneously, for ct > o-q , and the curve (/i')r makes 
its first appearance. From the asymptotic formulae (App. II) the latter 
runs to infinity between 7i and 0| , and now the solution curve which 
ended for Wi < ro < ji at Ji — {Ia)t is carried to infinity between I^ 
and (Ii)t . The asymptotic expression for /3^ versus a, given in formula 
(56) indicates that j3- tends to the isotropic value for the TMu mode. 
No further qualitative changes will take place in beha^dor of this mode 
as ro increases. 

As ro increases through U2 (the value at which the isotropic medium 
supports the TE12 mode), the {Oi)t curve makes its appearance, an 
event accompanied by the escape of the uppermost incipient solution 
curve (the one ending at (Ia)t — h) to infinity. The escape takes place 
in the same way as that of the TEn solution curve as ro passed through 
ui . The newly escaped curve, of course, represents the TEig-limit mode. 
The end points of the remaining incipient solution curves also jump 
discontinuously to their next higher neighbors as they did at ro = ih . 
The course of events as ro is increased further should now be abundantly 
clear, and is summarized in Table I on page 642. 

We now turn to the region < o- < tro and consider first the situation 
< To < Wi . It is clear that in the area bounded by 0„, Ou , (Ia')t and 
(7'B)r both G functions are negative. There is no simple geometrical argu- 
ment which determines the existence of a solution curve in this region. 
It is therefore necessary to use a type of analytic argument, which is use- 
ful in a number of other cases, although fully discussed only in the pre- 
sent instance. 

We show that the least value attained by G(Xi , a) in the admissible 
region for p = (which contains all regions admissible for other p-values) 
is greater than the maximum value of G(\2 , cr) in the range — 1 < X2 
< 0, ff > 0. Consider the variation of G(Xi , a) as the point Xi = 1, 
ff = 1 is approached along a line of constant x 'n the admissible region 
for 7> = (see Fig. 4). We have the relation 



G(\a) = \ 

X- 



l F(rox) - 1 
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For X negative, F(rox) is positive and thus, as X approaches unity from 
above G decreases. Again for x positive and ro < k.i , F{hx) is positive 
and thus, as A approaches unity from below, G also decreases. Thus for 
any x^ Gi^, a) takes on its least value at <t = 1, X = 1 and this value is 

\\nr,x) - 1]. 

The minimum value of this limit in this region is F(ri)) - 1 and is greater 
than —1. In the region - 1 < X2 < 0, a > 0, x^ is between and l,and 
the X curves run from £r = Otoff= '^■> . G will clearly decrease as <t 
increases from zero on any one of those curves. Thus G attains its maxi- 
mum on tr = 0, where its value is 



1 - V L I X I J' 



Since F{n-\/\ - V) is positive for ro < uj and | X | < 1 , G is clearly less 
than -1. In passing we note that for ro = Wi, both G functions may 
attain the value —1. 

As ro passes through Ui , the (Oo')r curve appears in the region under 
discussion and together with {Ib)t delimits the region carrying the TEn- 
solution cui-ve already discussed at length. No ciualitative changes occur 
iu that curve as ro is increased indefinitely. When ro exceeds ji , the 
(/i')r curve appears between (Oo')r and (Ia')t- Between (Ia')t and 
(Ii)t the G functions have a region of common sig n, ye t no solution 
curve arises there for a given p until ro reaches jVVl - p''* From then 
on, the 7i curve cuts (I/)t , see Fig. 9(c), and a solution curve exists 
between (/i')r and 7i . It is cut off at the intersection {I/)t — h; there, 

if = Q and 0-, — are given by the same parametric formulae (16-8) 

da 
applying to the cut-off of incipient modes, the parameter d being nega- 
tive. The curve begins at a = 0, where it satisfies the usual etiuation, 
which for this radius has two solutions. The solution \\-ith the smaller 
X, belonging to the present curve, tends to the isotropic TMu-limit as 
p __^ 0. At a fixed n , sufficiently below U2 , this mode does not exist at 



* There are some exceptions to this statement. When 4.82 < ro < H2 = 5.33 and 
V exceeds \/l - j'lV'-o", a double-valued 0= - ff curve exists between two positive 

a values For values of co still closer to v-. further regions of common sign may 
arise as a result of the interplay of the (Oi'JT and (Li')r curves. We hiive not 
examined these regions eloaely. Such duliious regions are confined to the immedi- 
ate neighborhoods below the w„ . 
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all when 



If fo is greater than ks , the (Oi')r curve has appeared. A new region of 
Uke signs of the G's arises between it and (/i')r , see Fig. 9(e), and con- 
tains a sohition curve. This ends at o-u , Xio and Ijegins at o- = at a value 
of Xi pertaining to the T^Mu-mode. Thus, it is clear that as i\, passed 
i(-2 , tlie end-point of the TMn curve jumped discontinuously from 
(/aOt — 1 1 to (To , Xio - This jump is anticipated as r^ approaches u^; the 
iS" — 0- curve first bulges beyond {Ia')t — h towards its later course and 
returns to that point with positive slope. As ?o increases further no change 
occurs in the qualitative behavior of the mode. It may be noted that 
above ih the mode exists for all p. 

Beyond Va = vz , at least part of the area between h and (Ia')t is an 
admissible region and does in fact contain the TE12 -solution curve. It 
begins at 0- = and Xi given by that solution of eqn. (39) which is, in 
the limit p = 0, the TE12 solution. It is cut off with ,3" = at {Ia)t 
- h , the end point relinquished by the TMn-solution curve. As rt, 
passes J2 , the TE12 solution retains its cut-off point, but, beyond n = Ws , 
it will transfer this point discontinuously to (To , Xio . Thereafter its course 
remains essentiaUy unaltered. Tables I, II and III show the progression 
of cut-off points of the various modes. 

It may be recalled that in the analysis of id < ro < Ji , the modes in 
17 < — ffo followed essentially the same course as in a > aa . This is also 
true of their progress with changing radius and of the escape process. 
The singular character of the (Oo)r curve and the presence of /« lead to 
some local changes in the progress of the modes but have no effect on 
their more salient features in this particular range of a. The scheme of 
progression of the end points is shoA\n in Table I. 

In contrast with the state of affairs in the region just discussed, the 
mode structure in the area between a = and a = — ao is very mark- 
edly affected by the presence of (Ou')t and ///. 

Wlien ro < «i , a solution curve exists l^etween cr = —1 — p, and 
a = -<Ta . It starts with /iJ" = at the intersection of (7^)7 and (V) with 
a slope given by (41). For sufficiently small ro , &^ tends to infinity as 
(T — * (Tu , smce the solution curve approaches the line 0/ or (O^)?- . Its 
shape is then given by (52), see Section 4.17. As r^ increases, Oo falls 
steadily. Eventually, for suflficiently large p, its minimum falls below 



• See footnote on page 628. 
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— o-Q . (0(,)t now has two branches for a > — <ro , which pass through the 

point (7 ^ —<To,M= —— making there a finite angle with each other. 

{Oa)T is completed by a loop in a- < — o-q , Fig. 9(g), which does not 
affect the incipient modes appreciably. The mode in question now has 
two branches. The first starts as before and ends at o- = — tro , where the 
associated ^^ is given by /3„- = Xq/o-o and X^ is the smaller root of 



It resumes at o- = — o-q and ^^ ^ ffi,^ ^ Xb/o-u , where Xa is the larger 
root of the above equation, progresses to smaller | <r |-values and then 
back to (T = — ffo where /3^ tends to infinity again in accordance with 
(52). Beyond ro = Wi , where Oo rises steadily from o- = — «= too- ^ with 
increasing X, one branch of (Oo)r m —ao < cr < disappears and only 
the second branch of the mode remains. Neither branch has an analogue, 
in ordinary waveguides; as p — > each lies in a smaller and smaller 
neighborhood of a- = 1, and finally vanishes into tr = 1, X2 = — 1. 

For To between i(i and ji there is a single solution curve starting at 
ff ^ and ending with ,S" = at {Ia)t — Is' ■ This may be identified 
in the limit p = 0, with the TEn-limit mode, and has already been fully 
discussed for Ui < rn < ji . No change in the formula for its cut-off point 
occurs up to To = W2 . A useful spot-point {Is — (Ii)t) along its course 
can be found when 



< 2> < 1 - 



iA^ 



and is given by (60). 

In the range ji < vq < «2 , a further solution curve (corresponding to 
the TMii-limit mode) can arise, provided 



V U2 



The radius at which it mil then first appear is 



'■0 = 



vr^ 



p- 



It begins on o- =0, according to (39) and is cut off, with fi^ =^ 0, at 

{Ia)t - h'. 

As To passes 112 , the cut-off point of the TEj solution curve moves dis- 
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a- 
Fig. 13 — Approximate course of modes in the ^- — w phuic for various p values 
audat two values of To . Fig. 13(a), above, r^ = 2.75, TEi, -limit modes; —1 <<7- < 1. 
Fig. 13(b), To = 2.75, TEu-Hmit mode and incipient TM,, -limit mode, <r > 1 . Fig. 
13(c) , To = 2.75, TEn-limit mode, o- < - 1 . Fig. 13(d) , (c) and (f) , m = 5.75, TE,,- 
limit and TEi2-limit modes; Fig. 13(d), -1 <ff < 1; Fig. 13(c), <r > 1; Fig. 13(f), 
ff < -1. Fig. 13(g), 13(h) and 13(i), r„ = 5.75, TM„-limit modes; Fig. 13(g); 
— 1 < ff < 1; Fig. 13(h), o- > 1; Fig. I3(i), <t < 1. It should be noted that a scale 

linear in -. — : — i is used for convenience when | a [ > o-q . 



1 



3^ 



continuously to o- = — tro , ^2 = — — and, simultaneously, the cut-off 

Co 

point of the TMn curve occupies the position relinquished by the for- 
mer. The TMn mode now exists for all p. A new solution curve (TE12) 
appears in the region })Ounded by 0/, (Oi)t and //, if p is not too large, 
terminatmg at {Ia)t — h' , the point left by the TMi terminus. (If p 
exceeds -\/l — ji^/-u/ this curve will not exist at all.) 

Figs. 13(a) to (i) and 14 show the approximate course of the ^ — a 
curves for the TEu mode at tq = 2.75 and for the TEu , TMn and TE12 
modes at r^ = 5.75. The incipient TMu mode at ro = 2.75 is shown 
for positive c, p only. They were computed by the methods outlined 
above. 

4.14. Guides of large radius. It is of some interest because of the high 
dielectric constant of ferrites to examine the behavior of the modes as 
the radius, tq , is allowed to become very large. The two sides of the 
G-equation will remain determinate for unHmited ro provided 
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Fig. 13(b) —See Fig. 13. 
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Fig. 13(rl) —See Fig. 13. 



remains finite, while 



1 - Xs^ 



or To 



1 - Ai^ 



1 — (xXn \ r 1 ~ ""Xiy 

becomes infinite imaginary. Examining the solutions obtained under 
these conditions it is possible to find expansions for /3 in inverse powers 
of ?■(! . These are as follows: 

for p > 0, a > I or p < 0, - ao < (T < 



/?'=! + 



P 



1 + 



p/2 



{43a) 



0- - 1 \ ' <T - ij n^' 

where the .r„ are the successive roots greater than zero of F{.Vu) = 1 
and the modes are associated with the x„ by the scheme: TEn — >■ Xi , 
TMn->.r2,TEi2^.r3,etc: 

for ;; > 0, < (7 < 0-0 or 7> < 0, 0- < -1 
V 



0- — 1 \ o- + 1/ ro- 



(43b) 
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with F(y„) = - 1 and TEn -^ y, ; TMn ~^ y^ ; TE12 -^ 1/3 , etc. 

for p > 0, (To < ff < 1 and p<0, —1<<t<— ao, 

there are no solutions. 

These formulae are valid for any p which is not itself so small as to be of 
order 1/n. If they are applied to modes varying as e^"*, where n = ±1 
and a, p are positive, they indicate the follo^'ing results: for cr > 1, 

n = ±1, j3" — *■ 1 + r; for 0-0 < a < 1, no n = ±1 modes; for 

P 



a - 1 



. These, in turn, may be classified 



in the follomug way. For o- > 1, ?i = — 1, and for < tr < co , ?t ^ +1 
which correspond to /i and k both positive, the propagation constant 
tends to the value for a plane wave whose direction of circular polariza- 
tion coincides with that of the wave guide pattern. For ag < a < 1, 
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where n is negative, n = dzl, no modes exist for large enough guide. 
For < a < ou , n = —1 and a > l,n ^ +1, the propagation constant 
tends to that for a plane wave whose polarization is in the opposite 
sense to that of the field pattern and here n is positive, but k is negative. 
An examination of the field pattern in this last case shows that most of 
the field energy is indeed associated mth a circular polarization opposite 
to that of the pattern as a whole. 

The discussion of the preceding sections shows that the complete 
structure of the mode spectrum for a guide filled \^ith lossless ferrite is 
very complex. It is also clear that for some combmations of guide size 
and magnetic parameters the course of an individual mode in the ^ — a 
plane may be quite involved. In particular, two values of ^' associated 
with the same mode often occur at a given a. The extent to ^\hich the 
complexity of the spectrum will be obseiTed in practice will depend 
principally upon the loss of the real ferrite and upon the guide radius. 
The effect of loss near tr = 1, where the incipient modes are crowded 
will be to cause simultaneous excitation of many of these and conse- 
quently a confused z dependence of the guide excitation. For values of 
ro just l^elow j„ , the point of escape of the TE modes, the latter exist 
over consideral)le ranges of a, see Fig. 12(a), and would probably be 
observable. The TE modes near /;„ also persist over a mde range, but 
are doulile-valucd. Concerning such double-valued waves it may be ob- 
served that from the results of the subsequent treatment of losses, it is 

clear that if ti — ; > 0, it is necessary to put the source of power at the 

opposite end of the guide. 

4.15. Losses, Faraday rotation and merit figure. So far the analysis has 
been ronccnicd with the loss-free mechiun. It is of some interest to 
determine the attenuation constant {the imaginary part of 0) that 
arises when losses are taken into account. As long as these are small, 
this can be done rather easily; in fact, sufficiently far from resonance 
{a = 1), for each formula giving /3", we can establish one giving the 
attenuation constant. 

If the lo.sses are of magnetic origin we utdize the fact (already demon- 
strated in section 2) that to first order in a, the pei-meabilities n, k are 
functions of tr + ja sgu p, and of no other combination of a, a. Since 
<T, a enter Maxwell's equations only through ^ and k, /3", which is derived 
from them, must likewise depend on <r through o- + Ja sgn p. Any formula 
for /3" derived for the loss-free medium can, therefore, be generalized to 
the lossy case by replacing u with c + ja sgn p, to first order in a. To 
this order, then, we find 
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(3^ = /3'^ + ja sgn p 



da 



where ^' is the propagation constant for the loss-free case. Thus 



jd - JP' - 



2^' a i cr I ' 



(44) 



and the last term on the right, (multipHed by our scaling variable 
00 = ti)-\/ju,eo) is the attenuation in nepers per meter. The present con- 
vention is that the waves propagate in the positive z direction, as exp 
(—j02). It follows that they will decrease in that direction only if 
a{j3')V^l ff I < 0. Occasionally this is not the case, and presumably 
indicates that the direction of the power flow^ opposes that of the phase 
velocity. 

For small dielectric loss, too, it is possible to derive formulae for the 
attenuation constant from those already obtained; obviously the latter 
depend on c only through e = eo — jei , and can therefore be expanded. 
But it must now be remembered that /3 was defined as ^„ctua\/^\/fi^e. 
and To as rnctuni coVTi^ so that the scaling parameter o}\/fi^e will make 
contributions to the imaginary part. It is then readily verified that 



2eo^'a(ro' 



iroY). 



(45) 



A few words may be said about the relation of Faraday rotation to 
the ^^ — cr curves. A linearly polarized plane wave traveling in the un- 
bounded medium along the magnetizing field can be regarded as the 




Fig. 14 — The course of the full}' developed modes (solid lines) and of some of 
the lower incipient modes (dotted lines) as a function of tr for ro = 5.75 and 
I p I = 0.6. 
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/ 

sum of right and left circular components which travel mth different 
propagation constants. If these are ^+ and ;S_ (measured in units of ^o) 
the plane of polarization of the resultant mil appear to rotate by 

X 1 

(/3+ — ,3-)/2 radians per reduced wavelength ^ = "^ ■ 

In the filled waveguide, on the other hand, it is no longer true that right 
and left circularly polarized modes add up to a plane polarized mode, 
as is readily seen by reference to the field components given in Appendix 
IV. To define Faraday rotation in a simple way it is therefore necessary 
to neglect changes in the field pattern due to the magnetization and 
consider only the changes in the propagation constants. Then the rota- 
tion of a mode wdth azimuthal mode number n will be ~ i^-i- — /3-). 

In the present case ?i = 1, and the 0+ , i3- are found from the curves of 
/3^ versus a, Figs. 13(a) to (i), for positive {a, p) and negative {a, p) re- 
spectively. 

The merit figure is defiiied as the ratio — radians rotation per neper 
loss — and is independent of path-length. For small losses, (neglectmg 
terms 0{a)), this ratio is 

Rl{d^ - /3-) ^1 ^+' - &J 

Im(l3+ + ,5_) ad0^ dp J ' 

a I I I '\ I I 
I cr I o I 0" I 

in the notation of the present Section. 

4.16. Formulae for the ferrite. 
1. Cut-off points 

Cut-off points win be classified into three types, 1, 2 and 3, accordmg 
to the nature of the intersecting curves which generate them. All points 
of a given type may be assigned an index which further identifies the 
generating curve. This will be written as a subscript. 
Type 1. Intersections of /„ - (Ia')t , a > 0, written as 1„ and of 
/./ — (/a) 7- , (T < 0, written as 1,/ 
0\= 
There is a parametric representation: 

|Xi,2| = e^ 



[ p I = 2 sinh d(l - ?; j and 

2 / 2\ 

II ^0 I / 1 '"0 ^ ^-s 



(46) 
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The slope at cut-off : 



d£ 

dtr /p 






3^ 



- 1 + 



^(ro)J 



■ sgn (T, 



e — e 



and 



a/3- 



p coth 6 \dcT Jp 

Type 2. Intersections of 

(0/)7- - 0^ at 0- = (To = 



(47) 



(48) 



^" + 1 - I ; type 2/ 



(0„)r — 0/ at ff -= —(TO ; type 2„ . 
Define Xo as a root of 



-1 



(To = 



^0 L 



1 - (1 - Xo') 



ro 



F-'0^o)J ]' 



(49) 



with the following convention: for 2„', Xo is negative and the 
n^ branch of ^^(Xo) is used; for 2^ , Xo is positive and the ri^ 
branch of F~^{\o) is used. Now 



/3^ = ^ 
Co 



Near cut-off 
'Xo 



To 



/3^ - '-^ + 

(To A 



1 + (70^ 



with 



(1 - Xo')A = 



2 Xo 



1 - 



Co 



Vq 



(To — (T 



1 — 0"Q^ y 1 + Xo ffo 



(50) 



(51) 



1 + Xo 1 (To 



((To -f- 2 I Xo I + (ToXo 



+ 



1 + (^o I Xc 



For 2o a special situation arises for 1 < ro < iii where there are two 
positive solutions for Xo . We write 2qi and 2o2 for the points corresponding 
to the smaller and greater of these. A special cut-oif point will be labeled 
2o«j and arises from 0/ and (Oo)r as Xi goes to infinity. For this point 
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Fig. 15 — Tho function 



Fir,) 



1 - F(ro 



related to cutoff of Type 3. 



7 ^ aa , p 



CO and near —ao we have 
I P I 



/3= = 



(2 - 1p 



0-0 



I)' 



(52) 



^(l"(| 0"|) I — I 

Type 3. Intersections of Ib — (J.i)t; for — o-q < o- < only; no sub- 
script is needed. 

/3' = 0, 



a = -1 -p, 
and 



1 i^(ro) 



(see Fig. 15) 



(54) 



(da)^ (dp), V 1 - F{n) 

The cut-ofT points of the modes follow various schemes in different 
ranges of a as indicated below. 

For a > (To we have Table I. When a < —an , 1,/ replaces 1„ in the 
Table I. "None" indicates that the mode exists, but has no cut-oft'. 

For < cr < CT() Ave have Table II. "N.P." in Table II indicates that 
the mode is not propagated. For — o-u < o- < we have Table III. In 
this range of o- one has alfio the mode without classical analogue. For 
n < Ml this is cut-otf at 3 and 2|ii and may have a second branch from 
2o2 to 2o« . For To > ui the second branch only exists. 
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Table I 



Mode 
Rndius 


TEii 


TMii 


TEii 


TNfn 


TEii etc. 


'"0 < Ui 


li 


1. 
\ 


h 


h 


U 


Ui < ro < ji 


None 


1. 


h 

I 


h 


U 


ji < ro < 1(2 


None 


None 


\ 


I3 


I4 


»2 < ro < j2 


None 


None 


None 


I2 


h 


J2 < ro < U3 


None 


None 


None 


None 


1. 


etc. 













II. Asymptotes, genesis of the modes and spot points 
For I a- 1 — J- a= there are asymptotic formulae : 
For TEin-modes 



P' 



1-^1 + 



ra' 



V 



For TMi„-modes 



f 






(55) 



(56) 



For ] 0- i > (To , all modes have their origin in the points o- = 1 , Xi = 1 
or a- = — 1,X2 = —1, where ^^^ oo. The variation of /SSvith tr in the 
neighborhood of these points is described by the two expansions: 

f or cr -^ 1 



+ 0(ff - 1) 



(57) 



where 

and 

with the scheme 



a„ + 1 ~ 2ro= 
F(x„) = 1 



a;„ > 0, 



ft 


1 


2 


3 


4 etc. 


Mode 


TEn 


TMn 


TEi, 


TMi2 etc. 
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Table II 



Mode 
Radius 


TEii 


TMu 


TEii 


TMu 


TEu etc. 


< 1-0 < III 


N.P. 


N.P. 


N.P. 


N.P. 


N.P. 


til < Co < Jl 


2o' 


N.P. 


N.P. 


N.P. 


N.P. 


jl < '■() < "2 


2/ 


1, 
\ 


N.P. 


N.P. 


N.P. 


"2 < ''0 < is 


2o' 


2i' 


ll 
i 


N.P. 


N.P. 


j- < rn < V:, 


2n' 


2/ 


ll 


U 


N.P. 


J/3 < ro < ji 


2,i' 


2,' 


2/ 


li 


h 


etc. 













Table III 



Mode 
Radius 


TEii 


TNrn 


TE,! 


TMii 


TEii etc. 


< ro < Ui 


N.P, 


N.P. 


N.P. 


N.P. 


N.P. 


111 < ro < j, 


3 

i 


N.P. 


N.P. 


N.P. 


N.P. 


>, < J-o < 11-= 


\ 


1/ 


N.P. 


N.P. 


N.P. 


n-i < - '■(! < J2 


2t 


3 


1,' 


N.P. 


N.P. 


ji < '■() < ^l■3 


2, 


3 

\ 


ll' 


I2' 


N.P. 


ii3 < rn < jl 


2i 


2.. 


3 


ll' 


I2' 


etc. 













for ij '^ — 1 
P 



^«- - 



1_ _^ J l - a„ _^ 2a£ /^ _ 2\ a^ 



where 



and 



!/" 



a„ + 1 2ro' ' 



+ 0(ff + 1) 



]• 



(58) 



nvn) = -1, 

with the same identification as above. 

For cr > (To a spot-point is given by 7„ — {Ia)T with 
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V 



^^ = Xi = 1 + 



1 + cr' 



and 



(59) 



1 ^ 2 - 2 

1 " Xl ^ 3n_ 

1 — o-Xi nr ' 



The identification scheme is again that shown above. 

For — (To < (T < an isolated identifiable point arises from Ib- 
{li^T which is expressible in inverse form by the relations 

3 



1 - 






+ ;^r and 



V = i0'- 1) 



1 -^ 

1 + ?^^^+ ^" 



U 



0' 



(60) 



for 



1-4 <^= < W^l_i!L 



The identification of the modes with n proceeds as in the earlier parts of 
this section. 

III. Small p. 

To order p^ there exist the following expansions for (f: for the TEi„- 
mode 

0' = &.' + 



1 - 0-2 



w„- - 1 



?> + 



2^0 



(1 ^ ^^)^ 



^0 



+ 



5 + 5m„ — 2w„ 



(^„^- 1)^ 4(U„^ - 1) 2 ..,^ 

(^0^ - 1) (1 - u.^) ^■■•' ^^'^ 



where 



For the TMi„-modes 



Q' = /3o^ 









1 - cr^ 



3,3o' - 2 2 



(1 - (7>')=2(1 - jSo^) "^ "■' 



(62) 
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where 

. 2 

The radius of convergence of these series is not known. It is clear that it 
wiU depend on a and will become smaller as a^ ^ 1. 

4.2. The Plasma (pn = 0, v„ = 1). The characteristic equation (26) 
may now be written 

-i-„ [XiF„(xi/-o) - »] = — „ [X2F„(x2ro) - n], (63) 

xr X2^ 

where (in contrast with the ferrite case) Xi,2 ^ /3Ai,2 . The X satisfy 

^^^2 _ (..-l)(l-gV..)-^.P.\^^ - /3^ = 0, (64) 

and the x's are given by 

Xi.2^ = (1 — /SVi'fl) — pffXi,2 (65) 

From the equations for pe and fE in terms of o-, q given in Section 2, 
equation (64) may be written 



Xi,.^ - ( n^ - ^^" ) ^u - /5' = 0, (66) 



or 

^2 



XiXs = -/5 (67a) 

(67b) 



Xi + X. - -j-^ -<r/3', 
1 - q' 



"" + (rXiXn . 



1 -5^ 

Elimination of if between equations (64) and (65) enables us to express 
Xi/ solely in terms of X1.2 , ph and pg • 

2 1 - Xl,2' 

X1.2 = - 



1 ^ - V"^ > 

I Al,2 



which, from the plasma formulae for ps , ve , can be written 



2 1 — X1.2 



x./ - f-^^ . (68) 

1 — (rXi,2 



/ 
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With these expressions for the x, the characteristic equation (63) takes 
the form 



H{\i , <T, n) = BO^-i , cr, To), 



where 



//CX,cr,ro)-^-^ [^^"('-"V^ 



- X^ 



- a\ 



— n 



(69) 



For given o-, and q, equations (67b) and (G9) are simultaneous equa- 
tions for Xi , X2 . When Xi , X2 have been found, fi"^ = -XiXi is known. 
Since 0^ must be positive, Xi , X5 must have opposite signs. As in the 
ferrite, the convention Xi > 0, Xa < will be adopted. Equation (67b) 
will hereafter be called the plasma relation. The transformation 

Xi -^ — X2 , X2 — > — Xi , a -^ —a 

leaves the plasma relation unchanged and changes n to -n in the 
H-equation. As more fully explained in connection with the ferrite sec- 
tion, it is therefore necessary to consider positive n only if a is allowed 
to take on negative as well as positive values. As before, only the first 
azimuthal mode number (n = ±1) is con-sidered in this paper. 

The method of analysis is the same as that used for the ftirrite. Here 
we shall only sketch the most important .steps; the reader Avill have no 
difficulty in completing the analysis Ijy referring to Section 4.11. For 
fixed ro , a contour map of H is drawn in the X, a plane (see Fig. 16 
draw^l for ro ^ 2.2). The gross features of this map arc determined 
by the fines H = 0, /f = ± 00 . For greater detail recourse is had to the 

lines — — = constant, along which values of H are readily generated. 

1 — ffX 
Further help is obtained from a knowledge of the location of the saddle 
point of H. The infinity curves are given liy the same formulae as for 
the ferrite, except that the line X ^ is no longer an infinity line: along 
X = 0, // = — 1. Zero curves are given by 

1 rijl - \') 
"" X \[F-\\-^)r 

The branches of ffX = 1 are also zero curves in the same restricted 
sense as for the G function. In the same notation as for the ferrite, all 
In curves pass through tr - 1, X = 1; all /„' curves through -], -1. 
The same is true for aU 0„ , 0„' curves (n > 0). The oidy exception is 
denoted by Oo , it arises from that branch of F~^ along which F~\i) = 0. 
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2.5 



2.0 - 



1.5 - 



1.0 - 



as - 



(T 



-2.0 - 



-2.5 




Fig. 16 ■ — The division of the X — o- pliiiie into regions of positive and negative 
H by the first few and / curves. Dotted regions are positive. Cross-hiitclicd re- 
gions contain tlie Iiigber O and I curves. 



Along all lines, a = \c -\- d, H tends to infinity except A\lien c = ra/u^ 

_ 2 

(the slopes of the linear asymptotes of curves). Along a ^ — :, \ -\- d, 

H tends to a value depending on d. More about the general behavior of 
// can be deri\'ed by means entirely analogous to those employed in the 
study of G. 
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X-pairs determined from the fi^-diagram do not necessarily solve the 
problem, since for a fixed q, they may not satisfy the plasma relation 
(G7b). To take it into account, we interpret it as a transformation of the 
whole of the second quadrant onto part of the first, and of the whole 
fourth quadrant on part of the third. Writing (67b) in the form 

2 
'^ -1 



<T tr 1 — 0-A2 

we see that the curves ^2 = const, transform into a bundle of hyperbola 
passing through the intersection of o-^ = 1 — g^ mth a = 1/X, that is, 
through 

1 . 

^" " a/T^^ ' "^o = V 1 - q'' 

These hyperbolae have vertical asymptotes Xi = — — , and cut the line 

A2 

ff = in — Xa . For a fixed positive <t < o-q , Xi decreases from l/c to a-/ 
(l — g^) as X2 increases from -co to 0, but, when cr > (Tq , Xi increases 
from l/tr to (r/(l — q^). Thus the second quadrant transforms into the 
region between a = X{1 — q^) and tr = 1/X in the first quadrant. Simi- 
larly, the inverse transformation X2 = T{\i) transforms the fourth 
quadrant into the region between o- = X(l — /) and a- — i/X in the 
third quadrant. Points outside these regions cannot be site of acceptable 
solutions of the H equation. In order to locate acceptable solutions, the 
H = equation is now written in the form 

H(\i , (T, ro) = i?(rCXi), <r, Tq) 

when 0- > 0, and in the form 

H{M,<T,n) = H{T(K2),<y,ro) 

when ff < 0. These equations represent the curves of intersection of 
the //-surfaces. Along each such curve, both /f-equation and plasma 
relation are satisfied. Their projections onto the first (or third) quadrant 
give Xi (or X2) as a function of tr, and hence X2 (or Xi) from the plasma re- 
lation. Thus )3^ = — X1X3 is kno^Ti along each solution curve. The rough 
locating of the solution curves, and the estabhshment of precise analyti- 
cal formulae near special points on them proceeds in complete analogy 
with the ferrite case. Here we shall consider only the radius ro '^ 2.2, 
as typical of radii large enough to permit propagation of the TEn mode 



1.5 



O.S 



GUIDED-WAVE PHOPAGATION THROUGH GYl{OiL\GNETIC MEDIA 64:9 

W, 





TE, 



l-Oi 



p'TE, 



./TE, 



-1.5 



-1.0 

A.3 



Fig. 17 — Geometrical exploration of solution curves for the plasma. The 
doited regions are e.xeluded by the plasma relation; the cross-hatched regions are 
those in whicli H{\, a) and HiT{\), a) have unlike sign. Solution curves may lie 
only in unshaded parts of the 1st and 3rd quadrants, g^ = 0.25, To ~ 2.0. 

through the unmagiietizecl plasma, but too small to admit higher modes. 
The solution curves for ?o '^ 2.2 are indicated roughly in Fig. 17. 

In the first quadrant, for a > ao , a. solution curve starts at Xi = 1, 
a = 1, passes through the intersection of Ii , (Ib')t and pi-oceeda to 
infinity as indicated. The formulae in Section 4.21 describe the cor- 
responding curve near cr = 1, and at o- — ?■ co, showing that the solution 
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curve describes the TEn-limit mode. Incipient modes also exist, just 
as in the ferrite; their end-points on a = (I — f/)Xi (or, briefly, on 
(Xa = 0)r) are now the points for which 7/(Xi , a) ^ —1 and, simulta- 
neously, ff = (1 — q)\i . 

Below (TO , there is only one solution curve for ra '^ 2.2. It begins at 
0- = 0, Xi = /3iao (= — Xa , by the plasma relation), where pi^o is the 
propagation constant of the TEu mode in the unmagnetized plasma. 
(In contrast with the ferrite, the plasma becomes isotropic as (t — * 0). 
It is cut off at the intersection of the contour //(Xi , o-) = — 1 in that 
region with (Xa = O),- . At that point |S^ = and cr ia best stated, thus: 



<r = (I - q-)V{l+ y')/[l + (1 - q')y% 
where y is the (unique) real root of 



njyro) = V{i 4- y')[i + (I - q')y']- 

Alternatively these two equations may (by varying y) be used to generate 
ro's and the corresponding cut-off values of (t. Of course, the two equa- 
tions are merely a re-statemcnt of the equations //(Xi , a) = —1,0- = 
Xi(l - q^), heed being paid to the fact that theargument of F is imaginary 
in the region considered for the radius under discussion. 

In the third quadrant for a < -ao, we also find the TEn-lknit mode. 
Its solution curve begi:is atX-^ —1, a = —1, and proceeds to o" = — » 
without passmg through any easily computed intersections of / curves. 
Formulae pertaining to the TEn mode in this range are stated in Section 
4.22. Again the incipient modes are found in their usual region. For 
> (T > — (To , the solution curve corresponding to the TEu mode 
begins at o- = 0, Xa - -0uo{= -Xi) and is cut off at the intersection 
of H(X2, 0-) = -1 withX2(l - g') = ff (or (Xi = 0)t). At that point 
/3^ = 0, and a is given by 



c= -a-Ql V(l - f-)/[l - (1 - q')y% 

where y is the least real root of 

F(r,y) = -V(l ~y')[l - (1 -q')y']- 

Alternatively, this equation can be used to generate ro , and the asso- 
ciated ff, if y is regarded as a parameter, ^vhich for ui < ro < ji is betw^een 
zero and unity. 

At a fixed ?-o the higher roots of the last equation with sign reversed and 
the corresponding a are associated with the cut-off's of the incipient 
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FiK- 18 ~ Cutoff value of a for the TEn-Hniit mode in the plasma as a function 
of Vti lor various q-. 

modes in the lower half plane. Similarly the real roots of 



f{m) = - V(i - y')[i - (1 - q')y'], 

and the corresponding 



(T == +(1 - ?) V(i - y')/[i - (1 - q')yn, 

are associated with the cut-offs of the incipient modes in the upper half- 
plane. These equations have been solved for the TEu mode and their 
solutions shown iu Figs. 18(a) and 18(b). 

4.21. Some formulae relating to the plasma {chiefly for TE-modes). 

The formulas given here employ dimensionless variables (Section 3) 
except where otherwise stated. 

Approximations for extreme values of a or q 
(a) small ff, q not near unity 

TEi„ mode: 



&' = pJ'VA^a-\-By + 



(70) 
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where 

3^ = 1-}^ 



Bm = — 



2 r .2 



(1 - q'yi uj - 1 \\2u 



i + .^^U^.-i)^/ 



"^ w- 1)^*^ 



2uJI]\ ' 



TM modes show no first order variation with a. 

(b) small q^j a not near unity. Here ^a and r^ are the actual propaga- 
tion constant and radius, without scahng factors : 
TEim mode : 



3 

/3„ = (1 - gjtujuoeo - —z- 



+ --ri.C-:;^,-0 + otf). 



(71) 



TMi;„ modes: 



.2 2 1/ ■ 2 \ 

^a = (1 - g )w ^oeo - -^ - w eoMo , . I 1 - V^ s 1 ■ t'2) 

ro^ 1 — ff2 Y coVoforoV 

(c) Approximation for large ff ; q not near unity. 
TEi^mode: 

R' = ^ - 1^' - ^g' ^ r?^^ 

'^ 1 - g^ ro^ (r(l - g^) (u^' - 1)" ^ '' 

All formulae in (a), (b), (c) apply to both positive and negative a- (right 
and left circular waves). Formulae 70 and 73 show that the first order 
changes in /3^, whether due to very large or very small tr, have coefficients 
that differ only in sign. 

The TEn mode near resonance (g < 1) 

Near a = +1, 

q' W " 7 (74) 



^=- 



1 - q'' 0- - 1 
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Near o- = — 1, 

^'-- '^' ^ (75) 



Cut-off of the TEn mode {q < 1) (i^i < n < ji) 
a positive: 

^' = 0; <7 = (1 - g') V(l +ymi+ (1 ~q')y'h (76) 
where y is the only real root or the smallest imaginary root of 

FUm) = V{i + y')[i + (1 - q')y']- 

a negative: 

H' = 0;<r= -(1 - q') V(l ~ y')/[l - (1 ~ qWh (77) 
where y is the smallest real root of 

Firoy) = -V(l - y')[l - (1 - q')y']. 
(See section 4.2 for further explanation). 

APPENDIX I. THE f -FUNCTION 

The function F{x) has been defined by the equation 

Using the infinite product for Ji{x) and differentiating logarithmically 
one finds 

Fix) = 1-21; ^T^. . (78) 

where Jiijn) = 0. Near one of its poles, j„ , F(x) behaves asjn/{x ~ jn). 
It is also useful to know the form of F{x) near one of its zeros, ?(„ , 
which are also zeros of J\{x). Such an expansion may conveniently be 
found by using the Ricatti equation satisfied by F{x), which is 

a:^= 1-.-^^-;^^ (79) 

ax 

The expansion near «„ is then 

F{u. -^y) =y\l; - ''"] -\ ^ + l] + ^^Sher terms. (80) 
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The equation may also be used to furnish an expansion near x = 0. 
This is 

2 4 

Fix) ^ I -~^^ + higher terms. (81) 

4 90 

Finally, putting x = jy, one finds from Eq. (79) for large y 

1 3 1 

F(jy) =y— --l----f- higher terms. (82) 

2 o y 

APPENDIX II. INFORMATION PERTAINING TO THE CONSTRUCTION OF 
G^-DIAGRAMS 

The accurate construction of the contours G = const, is conveniently 
based on the contours (l - ^^)/(l - o-X) = const, along any one of which 
G is a function of X alone. These contours are shown in Fig. 4. Their 
asymptotic properties are almo.st self-evident. 

The curves G ^ g = const, have various asymptotes. These, together 
with their range of \'ahdity, and their Folder transforms where needed 
are stated in Table IV. 

The formulas given in Table IV show that the curves G = const, 
generally have two kuids of asymptotes; linear and hyperbohc. Formula 
(83) shows the behavior of G along a line of constant finite slope unequal 
to n^/jn, the asymptotic slope of the /„ curves. Parallel to a line of 
slope n/jn oil G contours must be found, not just the restricted range 
given by the first formula. Writing a = ( r^/j^) X + a; in the equation 
G = g, and expanding F near its pole j„ , we find -i; m terms of g and ob- 
tain (84) which holds for all g, from - co to + o) . When j? = it also 
gives the linear asymptotes of On , 0„' curves except Oo , as is readily 
verified from the equation 

_ 1 T,\\ - \') 

" X x[^-Hx)? ' 

for the zero curves. 

Formula (85) shows how the G-contours tend towards o-X = 1 from 
the side ctX > 1 as X ^ 0. 

Formula (86) relates the asymptotic behavior of the curves G = g 
to the zero curves, g = 0, for small X. All G curves approach zero curves 
arbitrarily closely as X — > 0. The only exceptions are the infinity curves 
whose form near X = is 



1 /. ro- 



^ = i{^-g) (B7) 
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When ro = Wn+i , the 0„ curve merges with the 0„' curve at 



-2 ro 



Similarly all of the contours G ^ g of the sheet to which On belongs 
merge with the corresponding G = g of the sheet to which 0^' belongs, at 

2(1 + g) ro' 



(88) 



(Wn+l^ - 1) Un+1^ ' 

These remarks apply to all 0„ , 0„' curves, Oo included. 
The Oo-curve, for largo X, behaves as 

1 - ro 



and so tends to o- from above for ro < 1 and from below when ro > 1. 
(In fact, for ro < 1, Oo lies wholly in the first quadrant; when Ui > ro > 
1, Oo cuts ff = once.) 

The saddle points of G are most easily found by considering G in the 
coordinate net formed by the curves x^ = const, and X — const. At a 
saddle point 



dG 
d: 
and simultaneously 



^ ^ \(\nr.x) -1 



ax 



= 



The only saddle points that might be missed in this way are points at 

which the two derivatives are not independent, that is points where the 

X contours have vertical tangents, and it is easily verified that no saddle 

points exist there. 

fiC 
Proceeding with the differentiations, we find that rr- = gives 

oX 

^(rox) = 



or 



and so — =0 gives 
5x 



rox = M« (89) 



? + '^ F'M = 
X X 
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or 

„ 2 _ 1 

(90) 









rtn» 


= ^ 
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corresponding o-„. 


are 


given 


by 












1 - 


Xna 




2 








1 - 


O'nsAns 




ro= 



(91) 

and are all positive. Thus all saddle points lie in the first quadrant. At a 
saddle point G = —tq/k^ and therefore it is the intersection of two 
contours G = —r^lu^ . Forn > 1, one of these obeys the asymptotic 
formula (83), the other is asj'mptotie to /„ and 7„_i (see Fig. 5), and 
obeys (84), with n and n — 1, near those curves. For n = 1, one of them 
still follows formula (83), but the two "arms" of the other are asymp- 
totic to 0- = 1/X and Ii , and so follow (85) and (84) with n = 1 re- 
spectively. 
Three further facts useful in the construction of G-diagrams are: 

Y \2 J, 2 « 2 

Alonff a curve :. r = "V, G equals — —2; thus the zero curves of F 

are contours of constant G. 
Along \ = +1, 



As ff, X — ^ 1 along (r — 1) = a (X — 1) ; 



G- "+1 



2 
As cr, X -^ — 1 along {a -\- I) = a(X + 1) ; 



(93) 



— "-^ W V«-ii)-] 



As for GiT(\), a, To) we have, in addition to the asymptotic formulas 
in the table: 

Ib' transforms into = 1 + — ; — - . 

0- + 1 

The intersection of (I„)t or (Z/)?- ^"ith ff = is given by 
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APPENDIX II. THE FIELD COMPONENTS 

The field components are given here for the ferrite and for the plasma. 
They are normalized in such a way that E^ takes a simple form. It should 
be noted that the X's appearmg in these equations are those defined ui 
Sees. 4.11 and 4.2 for the ferrite and plasma respectively and have a 
different significance in the two cases. 

We write 

Then, for the ferrite, 

E. = [Ai(r) - A-^irW' , 






^*-) l|^.ter)-^lV- 



r x-r 



\_ (KM _ A\.s 

r xi' I Xi "J r xi \ X2 J J ' 

n, = -,- ^ -L {f„(x.O - ^ (1 - x^) 



Xi 



-^ F„(x2r) - f (1 - X2O 
r xr I X2 

and, for the plasma, 

E. = [Ai(r) - A^ijW"' , 

S, = - i TM;) i, {(1 _ xi^)i^n(xir) + 7.X1} 
/3 L r Xi- 



^inO 



e'"' 
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E, = \^ \ {X, F„(xir) + n{l - xf)( 

r X2^ J 



and 



'^•^'■^ ^ {X.F„(xir) + n) - ^ ~ {X.i^»(x.r) + »}] e^"' , 



L ^ xr 



r xi- 



L r xr 



^=«i_|i-„feO + »x,l].'-"' 



X2' 
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